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SPATIAL JOINT SEARCHER AND CHANNEL 

ESTIMATORS 

BACKGROUND 



[0001] This application is related to the following United States Patent applications, all 

5 simultaneously filed herewith: United States Patent Application Serial 10/ , 

(attorney docket 2380-776), entitled "Multi-Dimensional Joint Searcher And Channel 

Estimators"; United States Patent Application Serial 10/ , (attorney docket 2380- 

796), entitled "Temporal Joint Searcher And Channel Estimators", United States Patent 

Application Serial 10/ _, (attorney docket 2380-797), entitled "Spatio-Temporal 

10 Joint Searcher And Channel Estimators", all of which are incorporated by reference 
herein. 

[0002] FIELD OF THE INVENTION 

[0003] The present invention pertains to wireless telecommunications, and particularly 
to apparatus and method for determining a channel estimate for use in reconstructing 
15 data symbols transmitted over a channel. 

[0004] RELATED ART AND OTHER CONSIDERATIONS 

[0005] A wireless telecommunications unit typically includes both a transmitter and 
receiver for communicating with other wireless telecommunications units over a 
communication link. For wireless communications, the communication link is typically 

20 over an air interface (e.g., radio frequency interface). As used herein, a "wireless 
telecommunications unit" with its "wireless telecommunications receiver" can be 
included in a network node (e.g., a radio access network node such as a base station 
node, also called Node-B) or a terminal. Such "terminals" include mobile terminals 
such as user equipment units (UEs), which have also been called mobile stations, and 

25 include by way of example mobile telephones ("cellular" telephones), laptops with 
mobile termination. Thus, terminals can be, for example, portable, pocket, hand-held, 
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computer-included, or car-mounted mobile devices which communicate voice and/or 
data with radio access network. Alternatively, the terminals can be fixed wireless 
devices, e.g., fixed cellular devices/terminals which are part of a wireless local loop or 
the like. 

5 [0006] As shown simply in Fig. 32, a wireless telecommunications system includes a 
transmitting antenna 2300T and a receiving antenna 2300R. Channel 2302 describes 
the relation between the transmitting antenna 2300T and the receiving antenna 2300R, 
including the wireless interface. A signal, typically modulated into pulses, is 
transmitted over channel 2302 from transmitting antenna 2300T to receiving antenna 

10 2300R. The signal can comprise a "symbol" or a string of series of symbols, depicted 
as "m" in Fig. 32. The signal can carry user data and/or certain control data (e.g., a 
pilot bit or pilot sequence). The signal m as transmitted by the transmitting antenna 
2300T is convoluted with a channel impulse response h of the channel, so that the 
received signal at the receiving antenna 2300R is m * h (e.g., m convoluted with h). 

15 The received signal m * his applied to base band processing functionality 2304 of the 
receiver where the received signal undergoes radio frequency processing. The data 
portions of the received signal are applied to a detector 2306, which may be, for 
example, a demodulator such as a RAKE receiver. 

[0007] Most modern detectors attempt to recover a symbol estimate m from the 
20 received signal m * h . To do so, most sophisticated detectors expect to receive a 
"channel estimate" for use in modeling the channel over which the signal was 
transmitted. The accuracy of this channel estimate influences the accuracy and 
performance of the detector in estimating the actual symbol received over the channel. 

[0008] The modeling of the channel (which is necessary for most detectors) is 
25 facilitated by the control data, often in the form of a pilot bit or pilot sequence, which is 
transmitted by the transmitter. The control data, hereinafter referenced as "pilot data" 
for simplicity, is of a known or recognizable format or pattern. The pilot data is 
typically transmitted periodically by the transmitter source, and thus receipt of 
repetitions of the pilot data can be expected at the receiver at successive intervals. In 
30 view of factors such as relative motion of the transmitter and receiver, the successive 
intervals are not necessarily constant. The pilot data can be transmitted simultaneously 
with, or otherwise interspersed with, the user data. 
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[0009] In order to utilize the pilot data, wireless receivers typically include both a 
searcher and a channel estimator, such as searcher 2308 and channel estimator 2310 
shown in Fig. 32. For control data, the received signal m * his applied to the searcher 
2308, which determines a time of arrival (TO A). The time of arrival is then applied to 
5 the channel estimator 2310, which uses the time of arrival to determine the channel 
estimate h and then provides the channel estimate h to the detector 2306. Using the 
channel estimate h , the detector develops its estimate of the symbol, e.g., m . 

[00010] The receiver may receive the an original signal (e.g., short pulse signal) from 
the transmitter source through open space over a single, direct propagation path. 
10 Alternatively, in another environment having obstacles or other surfaces, the receiver 
may receive the same original signal over multiple propagation paths. In the multiple 
path case, the received signal appears at the receiver as a stream of pulses, each pulse 
having a different time delay in view of the corresponding propagation multipath over 
which the signal travelled, as well as possibly different amplitude and phase. 

15 [0001 1] Multipaths are created in a mobile radio channel by reflection of the signal 

from obstacles in the environment such as buildings, trees, cars, people, etc. Moreover, 
the mobile radio channel is dynamic in the sense that it is time varying because of 
relative motion affecting structures that create the multipaths, or due to movement of 
structures and objects in the surroundings (even if the transmitter and receiver are 

20 fixed). For a signal transmitted over a time varying multipath channel, the received 
corresponding multiple paths vary in time, location, attenuation, and phase. 

[00012] Some wireless telecommunications receivers capitalize upon the existence of 
the multipaths in order to achieve various advantages. Such receivers typically operate 
on the baseband signal to search for and identify the strongest multipaths along with 

25 their corresponding time delays. The receiver has a filter which operates on a power 
delay profile of the signal. The power delay profile can be conceptualized as a time- 
averaged refinement or other derivation of the channel impulse response. The searcher 
attempts to locate peaks in the power delay profile, each peak corresponding to arrival 
of a wavefront of the signal from a respective multipath. In many searchers the peaks 

30 also correspond to a channel tap of the filter. 
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[00013] A channel estimate h as applied to the detector therefore comprises a set of 
both time of arrivals (TO A) and complex channel coefficients, each pair of TOA and 
channel coefficients being associated with one of the arriving wavefronts. In other 
words, each arriving wavefront has a pair of members in the set, e.g., a TOA and a 
5 channel coefficient. The channel coefficients thus actually form a channel impulse 
response vector, so that the terms "channel coefficient" and "channel coefficients" as 
used hereinafter should be understood to refer to a channel impulse response vector. If 
there is only one wavefront, there is only one TOA and one channel coefficient in the 
set (one channel coefficient in the channel impulse response vector) . But for plural 
10 arriving wavefronts, there are a corresponding plurality of TOAs and channel 

coefficients. Ideally, the channel estimate h should provide as good an estimate of the 
channel impulse response as possible, thereby increasing performance of the detector as 
the detector makes its estimate m of the transmitted symbol m. 

[00014] The channel estimate is then supplied to the detector, such as RAKE type of 
15 demodulators. A RAKE demodulator typically allocates a number of parallel 

demodulators (called RAKE fingers) to the strongest multipath components of the 
received multipath signal as determined by the multipath search processor. In a 
wideband code division multiple access (WCDMA) radio access network, the outputs 
of each of the RAKE fingers are diversity-combined after corresponding delay 
20 compensation to generate a "best" demodulated signal that considerably improves the 
quality and reliability of the radio communications system. 

[00015] Conventionally, wireless telecommunications receivers first use their searchers 
to ascertain time of arrival of a wavefront. Subsequently, after the time of arrival has 
been determined by the searcher, the channel estimator utilizes the time of arrival to 
25 calculate a channel coefficient, which expresses both amplitude and phase of the signal. 

[00016] Some wireless telecommunications units have more than one antenna for 
receiving a same signal. In the prior art, the searcher attempts to locate peaks in the 
power delay profile for each antenna separately. In other words, for each antenna the 
searcher works more or less independently. See, for example, US Patent Publication 
30 US 2002/0048306, which is incorporated herein by reference. As such, the prior art 
searchers are essentially one dimensional. 
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[00017] As indicated above, the performance of a wireless receiver is considerably 
dependent upon the accuracy of the peak determination, i.e., time of arrival 
determination, performed by the searcher. The better the peak determination of the 
searcher, the better will be the overall performance of the receiver (e.g., less error rate). 

5 But in many instances it may be difficult for a searcher to find an actual peak in a 
power delay profile. As mentioned previously, in many searcher algorithms the peak 
corresponds to a channel tap. With such difficulty there is considerable risk of 
incorrectly choosing a peak. Moreover, it can then be difficult to estimate the actual 
channel tap value. Channels with low signal to noise ratios (SINRS) are particularly 

10 susceptible to these difficulties. 

[00018] What is needed, therefore, and an object of the present invention, is provision 
of apparatus and method for providing an improved channel estimate for a wireless 
telecommunications receiver. 

BRIEF SUMMARY 

15 [00019] A wireless communication receiver comprises an antenna array and a joint 

searcher and channel estimator. Plural antenna elements of the array provide respective 
plural signals (indicative of one or more arriving wavefronts) to the joint searcher and 
channel estimator. The joint searcher and channel estimator essentially concurrently 
considers the plural signals provided by the plural antennas for determining both a time 

20 of arrival and channel coefficient for each wavefront. The time of arrival and the 
channel coefficient are essentially concurrently determined by the joint searcher and 
channel estimator. The joint searcher and channel estimator applies the channel 
coefficient and the time of arrival to a detector which provides, e.g., a symbol estimate. 

[00020] The wireless communication receiver can be either a mobile terminal or a 
25 network node (e.g., a radio access network node such as a base station node, also called 
Node-B). In example illustrated embodiments, the antenna array can comprise a 
uniform linear array of plural antennas, but the joint searcher and channel estimator can 
also work with other types of arrays. 

[00021] Since it contemporaneously processes the signals from plural antennas over a 
30 sampling window in order to determine both time of arrival and the channel 
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coefficient, the joint searcher and channel estimator is considered a two dimensional 
unit. A first dimension is with reference to a time index of the sampling window, i.e., a 
sampling window time index. A second dimension is a spatial dimension imparted by 
the spacing of the plural antennas of the array. This spatial dimension, which involves 

5 essentially simultaneous and concurrent processing of signals from the plural antennas 
for the array in order to determine the time of arrival and channel coefficient, bestows 
on the joint searcher and channel estimator the distinction of being a "spatial" joint 
searcher and channel estimator. Thus, in the multiantenna embodiments, the spatial 
joint searcher and channel estimator is matched in a spatial domain to a direction of 

10 arrival. 

[00022] The spatial joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
embodiment the joint searcher and channel estimator includes a non-parametric type 
correlator (e.g., a correlator which performs a Fast Fourier Transform (FFT) 
15 calculation). In another example, illustrative embodiment the joint searcher and 
channel estimator utilizes a parametric approach. 

[00023] Concurrently using signals from all antenna elements of the antenna array, the 
joint searcher and channel estimator looks for pilot data in a sampling window, and 
generates the time of arrival and channel coefficient for each wavefront having a peak 

20 as seen in the sampling window. In so doing, for each sampling window the joint 
searcher and channel estimator stores the convoluted signals from the antennas. One 
example way of storing the signals for a sampling window is in a matrix, e.g., an 
antenna signal matrix. In constructing the antenna signal matrix, each of the plural 
antennas in the antenna array is represented by an antenna index. The joint searcher 

25 and channel estimator stores in the antenna signal matrix a complex value indicative of 
the signal received in the sampling window. The position or location of the complex 
value indicative of the signal received is determined by two indexes. The first index, 
conceptualized as being along the X axis of the antenna signal matrix, is the sampling 
window time index. The sampling window time index points to a time in the sampling 

30 window relative to a start of the sampling window. The second index, conceptualized 
as being along the Y axis of the antenna signal matrix, is the antenna index. 
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[00024] In an embodiment in which the joint searcher and channel estimator includes a 
correlator which performs a Fast Fourier Transform (FFT) calculation, the correlator 
considers a dimensional receptivity vector (formed from the antenna signal matrix 
using signals from the plural antennas of the antenna array for a particular sampling 
5 window time instance). The phase rotation speed, or frequency, of the dimensional 
receptivity vector for the sampling window time instance can be interpreted as the 
direction of arrival (DO A). There are plural possible frequencies for the dimensional 
receptivity vector, each of the plural possible frequencies corresponding to a possible 
direction of arrival (DOA) of a wavefront. The plural possible frequencies are 
10 represented by a frequency index. 

[00025] In conjunction with the Fast Fourier Transform (FFT) calculation, the 
correlator calculates Y(n,t) = FFT(n,X(:,t)), wherein t is the sampling window time 
index; X(:,t) is the complex antenna matrix (with the colon ":" representing all antenna 
indexes for one sampling window time index) and n is the frequency index. For a 
15 CDMA receiver, the correlator calculates Y(n,t) = Z Cj *FFT(n,X(:,t)), j = 1,K, wherein 
Cj is a coding sequence symbol value j and K is a length of the coding sequence. 

[00026] In the Fast Fourier Transform (FFT) calculation, the correlator output 
comprises Y(n,t). An analyzer of the joint searcher and channel estimator determines a 
maximum absolute value |Y(n,t)| max from the correlator output Y(n,t). A sampling 

20 window time index t_max at which |Y(n,t)| max occurs is chosen as the time of arrival of 
the arriving wavefront; a frequency index n_max at which |Y(n,t)| max occurs is chosen 
as the direction of arrival (DOA) of the arriving wavefront. An amplitude for the 
arriving wavefront is chosen by dividing |Y(n,t)| max by the number of antennas 
comprising the antenna array. A channel impulse response vector is determined from 

25 the directions of arrival for the arriving wavefronts. 

[00027] In another embodiment the joint searcher and channel estimator includes a 
parametric estimator which generates a parametric estimation output vector which is 
utilized by a channel estimate generator to generate the time of arrival and the channel 
coefficient. The parametric output estimation vector has a sampling window time index 
30 and a spatial parameter for each time index. The spatial parameter includes spatial 

frequency and spatial amplitude. The channel estimate generator uses spatial amplitude 
values of elements of the parametric estimation output vector to determine the time of 
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arrival and spatial frequency values to determine the direction of arrival of the arriving 
wavefront. Wavefronts in the sampling window are associated with each element of the 
parametric estimation output vector which has the sufficiently high spatial amplitude 
parameter value. The channel estimate generator uses a sampling window time index 
5 for an element of the parametric estimation output vector having a sufficiently high 
absolute value as the time of arrival of the corresponding arriving wavefront. The 
direction of arrival of the arriving wavefront is the spatial frequency parameter value of 
the identified time of arrival. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 [00028] The foregoing and other objects, features, and advantages of the invention will 
be apparent from the following more particular description of preferred embodiments as 
illustrated in the accompanying drawings in which reference characters refer to the 
same parts throughout the various views. The drawings are not necessarily to scale, 
emphasis instead being placed upon illustrating the principles of the invention. 

15 [00029] Fig. 1 is a schematic view of an example, generic wireless telecommunications 
receiver which includes a joint searcher and channel estimator. 

[00030] Fig. 2A and Fig. 2B are schematic views of differing example embodiments of 
spatial joint searcher and channel estimators, each shown with an antenna array. 

[00031] Fig. 3 is a diagrammatic view illustrating a signal emanating from a 
20 transmitting antenna along three separate multipaths to an antenna array of a wireless 
telecommunications receiver. 

[00032] Fig. 4 is a diagrammatic view of a wavefront travelling toward an antenna 
array. 

[00033] Fig. 5A and Fig. 5B are diagrammatic views depicting signals obtained upon 
25 arrival of a wavefront at an antenna array. 



[00034] Fig. 6 is a diagrammatic view of an antenna signal matrix. 
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[00035] Fig. 7 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatial 
joint searcher and channel estimator, with the matrix analyzer using a non-parametric 
analysis technique. 

5 [00036] Fig. 8A, Fig. 8B, Fig. 8C(1), Fig. 8C(2), and Fig. 8C(3) and are diagrammatic 
views depicting results of a comparative operational evaluation contrasting 
performance of a spatial joint searcher and channel estimator with a conventional 
searcher. 

[00037] Fig. 9A is a diagrammatic view of an antenna signal matrix; an antenna weight 
10 vector; and a non-parametric output estimation vector. 

[00038] Fig. 9B is a diagrammatic view of an antenna signal matrix and a parametric 
output estimation vector. 

[00039] Fig. 10 is a flowcharting showing representative basic steps performed by 
matrix analyzer and channel estimate generator of an example embodiment of a spatial 
15 joint searcher and channel estimator, with the matrix analyzer using a parametric 
analysis technique. 

[00040] Fig. 1 1 is a diagrammatic view illustrating a coherent combination of signal 
outputs by a joint searcher and channel estimator. 

[00041] Fig. 12 is a schematic view for illustrating how an antenna weight vector 
20 facilitates the coherent combination illustrated in Fig. 11. 

[00042] Fig. 13A is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
searcher and channel estimator comprising a matrix analyzer which employs a non- 
parametric analysis technique. 

25 [00043] Fig. 13B is a schematic view of an example embodiment of a temporal joint 
searcher and channel estimator shown with an antenna array, the temporal joint 
searcher and channel estimator comprising a matrix analyzer which employs a 
parametric analysis technique. 
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[00044] Fig. 14 is a diagrammatic view depicting a sequence of sets of pilot data and 
user data received by a receiver which utilizes a temporal joint searcher and channel 
estimator, as well as an antenna signal matrix utilized by the temporal joint searcher 
and channel estimator. 

5 [00045] Fig. 15 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a 
temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
parametric analysis technique. 

[00046] Fig. 16A is a diagrammatic view of an antenna signal matrix; a doppler weight 
10 vector; and a non-parametric estimation output vector for a temporal joint searcher and 
channel estimator. 

[00047] Fig. 16B is a diagrammatic view of an antenna signal matrix and a parametric 
estimation output vector for a temporal joint searcher and channel estimator. 

[00048] Fig. 17 is a flowcharting showing representative basic steps performed by a 
15 matrix analyzer and channel estimate generator of an example embodiment of a 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
parametric analysis technique. 

[00049] Fig. 18A is a schematic view of an example embodiment of a spatio-temporal 
joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
20 joint searcher and channel estimator comprising a matrix analyzer which employs a 
non-parametric analysis technique. 

[00050] Fig. 18B is a schematic view of an example embodiment of a spatio-temporal 
joint searcher and channel estimator shown with an antenna array, the spatio-temporal 
joint searcher and channel estimator comprising a matrix analyzer which employs a 
25 parametric analysis technique. 

[00051] Fig. 19 is a diagrammatic view depicting a sequence of sets of pilot data and 
user data received by a receiver which utilizes a combined spatial/temporal joint 
searcher and channel estimator, as well as an antenna signal matrix utilized thereby. 
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[00052] Fig. 20 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a non- 
parametric analysis technique. 

5 [00053] Fig. 21 is a diagrammatic view of an antenna signal matrix; a doppler weight 
and antenna weight vector; and a non-parametric estimation output vector for an 
example embodiment of a spatio-temporal joint searcher and channel estimator which 
operates in a three dimensional essentially concurrent mode. 

[00054] Fig. 22A and Fig. 22B are diagrammatic views depicting operation of a first 
10 alternative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 

[00055] Fig. 23 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 

[00056] Fig. 24A and Fig. 24 Bare diagrammatic views depicting operation of a second 
15 alternative implementation of a non-parametric, sequential spatio-temporal joint 
searcher and channel estimator. 

[00057] Fig. 25 describes the procedure of non-parmetric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

[00058] Fig. 26 is a diagrammatic view of an antenna signal matrix and a parametric 
20 estimation output vector for an example embodiment of a spatio-temporal joint searcher 
and channel estimator. 

[00059] Fig. 27 is a flowcharting showing representative basic steps performed by a 
matrix analyzer and channel estimate generator of an example embodiment of a spatio- 
temporal joint searcher and channel estimator, with the matrix analyzer using a 
25 parametric analysis technique. 

[00060] Fig. 28A and Fig. 28B are diagrammatic views depicting operation of a first 
alternative implementation of a parametric, sequential spatio-temporal joint searcher 
and channel estimator. 
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[00061] Fig. 29 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the spatial processing is followed by the temporal processing. 

[00062] Fig. 30A and Fig. 30B are diagrammatic views depicting operation of a second 
alternative implementation of a parametric, sequential spatio-temporal joint searcher 
5 and channel estimator. 

[00063] Fig. 31 describes the procedure of the parametric approach for spatio-temporal 
sequenced method where the temporal processing is followed by the spatial processing. 

[00064] Fig. 32 is a schematic view of a conventional wireless telecommunications 
receiver. 

io DETAILED DESCRIPTION OF THE DRAWINGS 

[00065] In the following description, for purposes of explanation and not limitation, 
specific details are set forth such as particular architectures, interfaces, techniques, etc. 
in order to provide a thorough understanding of the present invention. However, it will 
be apparent to those skilled in the art that the present invention may be practiced in 
15 other embodiments that depart from these specific details. In other instances, detailed 
descriptions of well-known devices, circuits, and methods are omitted so as not to 
obscure the description of the present invention with unnecessary detail. Moreover, 
individual function blocks are shown in some of the figures. 

[00066] Fig. 1 shows an example, generic wireless telecommunications receiver 20 
20 which, as mentioned before, can be included in a network node or a terminal, e.g. 
mobile terminal. The wireless telecommunications receiver 20 includes an antenna 
structure or array 22; a joint searcher and channel estimator 24; a detector 26; and, a 
timing and control unit 28. Optionally, as depicted by broken line, the receiver 20 may 
include a code sequence generator 30. 

25 [00067] As broadly employed herein, the antenna array 22 can comprise one or more 
antenna elements. Signal(s) from the antenna array 22 are applied both to joint searcher 
and channel estimator 24 and detector 26. The signal(s) from the antenna array 22 
comprise a channel impulse response vector if the antenna array 22 comprises more 
than one antenna element. 



[00068] In the likely event that the signal(s) have been encoded by, e.g., a spreading 
code or the like, both joint searcher and channel estimator 24 and detector 26 are 
connected to operate in conjunction with code sequence generator 30. The timing and 
control unit 28 generates timing (e.g., synchronization) and control signals which are 
5 provided to detector 26 and to joint searcher and channel estimator 24. 

[00069] It will be appreciated that the receiver may include, e.g., downstream from the 
antenna array, certain radio frequency processing functionality and radio frequency 
demodulating functionality, so that the signals applied to joint searcher and channel 
estimator 24 and detector 26 are baseband signals. The illustrated structure of wireless 
telecommunications receiver 20 of Fig. 1 thus essentially concerns processing of the 
baseband signal(s). 

[00070] Various non-limiting, representative examples of differing embodiments of 
joint searcher and channel estimators are described below. Ensuing descriptions of 
operation of wireless telecommunications receivers with these differing embodiments 
15 are premised on certain assumptions. Some of these assumptions are related to a 
channel model which conceptualizes electromagnetic fields as arriving in a discrete 
number of wavefronts at the wireless telecommunications receiver, and particularly 
arriving at one or more antenna elements which may be employed in the antenna array 
22. 

20 [00071] As used herein, a "sampling window" comprises consecutive time slots (or, in 
a CDMA system, for example, "chips") obtained from a given antenna and analyzed by 
a joint searcher and channel estimator. As described in more detail hereinafter, 
embodiments of joint searcher and channel estimators operate upon an antenna signal 
matrix formed from plural sampling windows. In some embodiments, hereinafter 

25 referenced as "spatial" joint searcher and channel estimators, the antenna signal matrix 
is formed from sampling windows obtained from plural antennas. In other 
embodiments, hereinafter referenced as "temporal" joint searcher and channel 
estimators, the antenna signal matrix is formed with respect to a single antenna, but 
formed from sampling windows obtained by that antenna for successive sets of pilot 

30 data (occurring over time). In yet other embodiments, hereinafter referenced as spatio- 
temporal joint searcher and channel estimators, the antenna signal matrix is formed 
both spatially and temporally. 
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[00072] For purposes of the technology described herein, the antenna array 22 is 
conceptualized as acquiring "dimensionally differentiated" signals. The joint searcher 
and channel estimator essentially concurrently uses the dimensionally differentiated 
signals provided by the antenna array for determining, for each arriving wavefront, both 
5 a time of arrival (TO A) and a channel coefficient. For the spatial joint searcher and 
channel estimator, wherein the antenna structure comprises an array of plural antennas 
having spaced apart or spatially separated antenna elements, the signals acquired by 
different antennas of the array are dimensionally differentiated with regard to a spatial 
dimension. For the temporal joint searcher and channel estimator, wherein the antenna 

10 structure comprises an antenna which provides signals for each of successive sets of 
pilot data received at separated time intervals, the signals acquired by the antenna are 
dimensionally differentiated with regard to a temporal or time dimension. For the 
spatio-temporal joint searcher and channel estimator, having both the antenna structure 
comprising an array of plural antennas and one or more antennas receiving the 

15 successive sets of pilot data, the signals acquired by the antenna are dimensionally 

differentiated with regard both to a spatial dimension and a temporal or time dimension. 

[00073] The joint searcher and channel estimators are, in some instances, said to 
perform a "concurrent" determination of time of arrival and some other quantity, e.g., 
direction of arrival or doppler shift frequency. In this sense "concurrent" means that the 
20 quantities or determinations could be derived in parallel from a result of an outcome- 
determinative operation, e.g. a non-parametric technique such as a Fast Fourier 
Transform or a parametric technique. 

[00074] SPATIAL JOINT SEARCHERS/ESTIMATORS 

[00075] In some embodiments, the joint searcher and channel estimator 
25 contemporaneously processes the signals from plural antennas over a sampling window 
in order to determine both time of arrival and the channel coefficient. In these 
embodiments the joint searcher and channel estimator is essentially a two dimensional 
unit, with a second dimension being a spatial dimension imparted by the spacing of the 
plural antennas of the array. This spatial dimension, which involves essentially 
30 simultaneous and concurrent processing of signals from the plural antennas for the 

array in order to determine the time of arrival and channel coefficient, bestows on these 
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embodiments of the joint searcher and channel estimator the distinction of being a 
"spatial" joint searcher and channel estimator. 

[00076] The spatial joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
5 embodiment the joint searcher and channel estimator includes a non-parametric type 
correlator (e.g., a correlator which performs a Fast Fourier Transform (FFT) 
calculation). In another example, illustrative embodiment the joint searcher and 
channel estimator utilizes a parametric approach 

[00077] Fig. 2A illustrates one example embodiment of a spatial joint searcher and 
10 channel estimator 24-2A which uses a non-parametric technique for determining time 

of arrival and channel estimate, as well as an associated example antenna array 22-2A. 

The antenna array 22-2A includes, by way of non-limiting example, four antenna 

elements 22-2A-1 through 22-2A-4. While the antenna elements 22-2 A- 1 through 22- 

2A-4 are shown as forming a uniform linear array (ULA), it should be understood that 
15 antenna configurations other than a uniform linear are possible, and that the number of 

antenna elements in the antenna array may vary (e.g., the number of antenna elements 

is not limited to four). 

[00078] There are coherency requirements for the antenna elements of antenna array 
22-2A, and for antenna elements for all other plural antenna arrays described herein. 
20 The coherency requirement can be fulfilled by the plural antenna elements being 

synchronized. Alternatively, even if the plural antenna elements are not synchronized, 
but their phase differences are known, the coherency requirement can be fulfilled by 
compensating for the known phase difference. 

[00079] The complex baseband signals obtained from the antenna elements are each 
25 applied to joint searcher and channel estimator 24-2A, as well as to a detector (not 
illustrated in Fig. 2A). The joint searcher and channel estimator 24-2 A comprises an 
antenna signal matrix handling unit 40-2 A. In one particular example manifestation, 
antenna signal matrix handling unit 40-2A includes antenna signal matrix generator 42- 
2A and antenna signal matrix memory 44-2A. A matrix analyzer, which for the non- 
30 parametric technique of Fig. 2A can be a correlator 50-2A, operates on complex values 
stored in antenna signal matrix memory 44-2A. The correlator 50-2A preferably 
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comprises a filter. The correlator 50-2A generates certain output values, which may be 
stored, e.g., in correlator output value memory 52-2A. The joint searcher and channel 
estimator 24-2 A further comprises a channel estimate (CE) generator 60-2 A. In the 
illustrated example embodiment, the channel estimate (CE) generator 60-2A comprises 
5 a correlator output analyzer 62-2A and a detector interface 64-2A. The detector 
interface 64-2A generates, for each wavefront, both a time of arrival (TOA) and a 
channel coefficient (CC). In Fig. 2A, the time of arrival and channel coefficient output 
by detector interface 64 are applied to the detector on lines 66-2A and 68-2A, 
respectively. 

[00080] In Fig. 2A, and other embodiments described herein, the transmitted 
electromagnetic signal is assumed to arrive at the receiver in a number of discrete 
electromagnetic wavefronts. A number of discrete electromagnetic wavefronts is 
presumed in order to accommodate the multipath phenomena discussed above. For 
example, Fig. 3 illustrates a signal emanating from a transmitting antenna 70 along 
three separate multipaths P u P 2 , and P 3 to antenna array 22. Each multipath has its 
individual amplitude, and accordingly has an associated complex number "a" of the 
baseband signal and a time delay x. For example, multipath Pi has associated complex 
number ai and associated time delay x\\ multipath P 2 has associated complex number a 2 
and associated time delay x 2 ; and so forth. As illustrated in Fig. 3, multipath P { is a 
relatively direct path between transmitting antenna 70 and antenna array 22; while 
multipath P 2 and multipath P 3 are reflected off obstacles 72 2 and 72 3 , respectively. 
Thus, the time delay x x for multipath Pi is shorter than the time delay x 2 for multipath 
P 2 , which in turn is shorter than the time delay x 3 for multipath P 3 . Similarly, barring 
other phenomena, it would be expected that the complex number a t for multipath P x is 
greater than the complex number a 2 for multipath P 2 , and so forth. 

[00081] For sake of discussion, the electromagnetic wavefronts are assumed to be plane 
("planar") electromagnetic wave fronts, such as the single wavefront 76 illustrated in 
Fig. 4 as traveling toward the antenna array. In all embodiments described herein it 
should be understood that the wavefronts need not be planar wavefronts, but that any 
30 other known form of wavefront may be considered in similar manner. Moreover, it 
should be kept in mind that Fig. 4 represents arrival of only one wavefront, but 
typically plural wavefronts are incident on an antenna array. 
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[00082] As further shown by Fig. 4, due to incidence of an individual wavefront, the 
output (e.g., signal) from each antenna element has its version of the complex number 
for the wavefront. For example, for a wavefront for the first multipath Pi of Fig. 3, the 
antenna element 22-1 outputs a complex number a M , antenna element 22-2 outputs a 
5 complex number a { _ 29 and so forth. The numbers are complex, and in the particular case 
that (1) the antenna elements are identical; (2) there is coherency, and (3) the plane 
wave has constant amplitude within the width of array, the absolute values of the 
numbers are the same. Furthermore, with respect to the same arriving wavefront, each 
antenna detects the arriving signal as having a phase. For example, for the wavefront 
10 for the first multipath Pi of Fig. 3, the output of antenna element 22-1 has a phase 8 M , 
the output of antenna element 22-2 has a phase 9i_ 2 , and so forth. 

[00083] Signals obtained upon arrival of a wavefront at a uniform linear array (ULA) 
antenna array are illustrated both in Fig. 5 A and in Fig. 5B. Fig. 5 A particularly shows, 
for each of four antennas 22-1 through 22-4, plane wave propagation over the antenna 

15 elements for a fixed time (chip) index, and resulting respective output pulses 78 (e.g., 
output pulses 78i through 78 4 ). For each corresponding antenna, Fig. 5B shows the 
pulse as a complex number and with the argument,9, of the complex number. The 
argument (0) corresponds to the phase of the received signal. The rate at which the 0 
values change (e.g., the rate at which the phase rotates) over time is known as the phase 

20 rotation speed, or frequency. The phase rotation for the wavefront with this array of 
antennas is depicted by the increasing angular value of 9 through the range of 9^ 82, 9 3j 
9 4 , and thus the frequency is the rate of change of this angular value over time. The 
phase rotation speed is constant. The speed of the linear phase propagation is 
dependent on the direction of arrival (DOA) of the incident wavefront. 

25 [00084] In the joint searcher and channel estimator 24-2A of Fig. 2A, the antenna 

matrix handling unit 40-2 A samples the complex baseband signals from each antenna 
element. Using the sampled complex baseband signals, antenna signal matrix generator 
42-2A generates an antenna signal matrix such as antenna signal matrix 80 illustrated in 
Fig. 6. The antenna signal matrix 80 may be stored in any convenient fashion, such as 

30 antenna matrix memory 44-2A. 

[00085] The antenna signal matrix 80 is a two dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 80 as a 
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function of two different indexes. For the antenna signal matrix 80 shown in Fig. 6, a 
first index is a sampling window time index, illustrated along the X axis of Fig. 6. For 
embodiments which utilize spreading codes or similar codes, the first index may be, for 
example, a chip index. Thus, the sampling window time index points to a time in the 
5 sampling window relative to a start of the sampling window. In the antenna signal 

matrix 80 of Fig. 6, a second index, shown along the Y axis, is an antenna index (which 
serves as a dimensional differentiation index). The antenna index points to a different 
row of the antenna signal matrix 80, each row being associated with a different antenna 
element in antenna array 22. Fig. 6 shows four rows in antenna signal matrix 80 for 
10 consistency with the previous examples of an antenna array comprising four antenna 
elements. It is reiterated, however, that the number of antennas in an antenna array, and 
thus the number of rows in antenna signal matrix 80 and the maximum value of the 
antenna index, can vary from receiver to receiver, and that the choice of four antenna is 
only illustrative for sake of example. 

15 [00086] The antenna signal matrix 80 is conceptualized as storing "dimensionally 
differentiated" signals acquired from the antenna array. For the spatial joint searcher 
and channel estimator, wherein the antenna structure comprises an array of plural 
antennas having spaced apart or spatially separated antenna elements, the signals 
acquired by different antennas of the array are dimensionally differentiated with regard 

20 to a spatial dimension. That is, for a given column of antenna signal matrix 80, the 
values in each row are dimensionally differentiated in the sense that they are acquired 
from different antenna elements which are separated in a spatial dimension in view of 
the separate physical placement of each antenna element with respect to other antenna 
elements of the array. 

25 [00087] For sake of simplicity, the complex values stored in antenna signal matrix 80, 
including the complex values obtained from the antennas, are not illustrated in Fig. 6. 
Such complex values would be illustrated in a third dimension, e.g., out of the plane of 
Fig. 6. The antenna signal matrix 80 includes both complex white noise and (for the 
sake of the present illustration) a complex sample for at least one wavefront (planar or 

30 other known shape). As stored in antenna signal matrix 80, the wavefronts have known 
phase (temporal, non-coherent detection), and are modulated code sequences. 

[00088] In conjunction with the antenna signal matrix 80 of Fig. 6, and particularly a 
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WCDMA case in which spacing of the antenna elements in the antenna array is not too 
far apart, the plane wavefront arriving at the antenna array can be considered to arrive 
in the same sampling window time index (or chip index). 

[00089] The complex values stored for each column of the antenna signal matrix 80 of 
5 Fig. 6 can be conceptualized as a dimensional receptivity vector. That is, a dimensional 
receptivity vector is formed with respect to a single sampling window time instance and 
with complex values from each of the plural antennas of the antenna array. Each 
element taken from a unique row of antenna signal matrix 80 has a different phase in 
the manner of the differing 0 values illustrated in Fig. 5. As received by the differing 

10 antenna elements, for the spatial joint searcher and channel estimator the change in 
phase over time is the frequency for the dimensional receptivity vector. If the wave 
arrives e.g. straight ahead, the angles could be the same. The phase rotation speed, or 
frequency, of the dimensional receptivity vector, for the sampling window time 
instance can be interpreted as the direction of arrival (DO A). Thus, each dimensional 

15 receptivity vector corresponds to a separate direction of arrival. There are plural 

possible frequencies for the dimensional receptivity vector, each of the plural possible 
frequencies corresponding to a different possible direction of arrival (DOA) of a 
wavefront. For the non-parametric techniques herein employed, the plural possible 
frequencies can be a continuous range of frequencies. For sake of differentiating the 

20 plural possible frequencies, the plural possible frequencies are each represented by a 
frequency index. 

[00090] The channel estimate generator 60-2A (see Fig. 2A) seeks to develop a 
"composite" channel estimate based on the complex values stored in antenna signal 
matrix 80. At this point it should be appreciated that, since the antenna array 22-2A has 

25 plural antenna elements, there are a corresponding plurality of channels through which 
wavefronts are received, and accordingly there could also be a separate channel impulse 
response or separate channel estimate for each of the plural channels. But by storing 
the complex samples in antenna signal matrix 80 in the manner aforedescribed, and by 
concurrently finding the time of arrival (TO A) and channel coefficients over the entire 

30 antenna signal matrix 80, the channel estimate generator 60-2A provides a channel 

estimate which encompasses all channels for all antenna elements and for this reason is 
known as a "composite" channel estimate. 



20 

[00091] The composite channel estimate comprises, as mentioned before, a time of 
arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 
window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
5 a time of arrival (TOA) and channel coefficient. The task for correlator 50-2A is thus 
to locate a value or "tone" in antenna signal matrix 80 that best corresponds to an 
arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 
sampling window. 

[00092] The task of locating a value or "tone" in antenna signal matrix 80 that best 
10 corresponds to an arriving wavefront can be accomplished by various techniques, 
including both parametric and non-parametric techniques. A Fast Fourier Transform 
(FFT) technique as discussed below is just one representative and illustrative example 
non-parametric type of correlator which can be utilized. 

[00093] Fig. 7 depicts example basic steps performed by an example correlator 50-2A 
15 and correlator output analyzer 62-2A in conjunction with the Fast Fourier Transform 
(FFT) calculation. As step 7-1, the correlator 50-2A of Fig. 2A calculates Expression 1. 

Y(n,t) = FFT(n,X(:,t)) Expression 1 

In Expression 1, t is the sampling window time index; X(:,t) is the complex antenna 
matrix (with the colon ":" representing all antenna indexes for one sampling window 
20 time index); and, n is the frequency index. Each FFT calculation is thus a one 

dimensional FFT calculation on the baseband signal, and corresponds to a specific 
direction of arrival (as depicted by the frequency index) and set of antenna weights 
which, in practice, are the FFT weights. 

[00094] The output of correlator 50-2A, i.e., the Y(n,t) values computed using 
25 Expression 1, are stored as correlator output values. The correlator output values can 
be stored, for example, in the correlator output value memory 52-2A of Fig. 2A. 

[00095] The correlator output analyzer 62-2 A of channel estimate (CE) generator 60- 
2A searches the correlator output values and (as step 7-2) determines therefrom a 
maximum absolute value |Y(n,t)| max . This maximum absolute value |Y(n,t)| max is 
30 utilized by correlator output analyzer 62-2A to determine both the direction of arrival 
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(DOA) and time of arrival (TOA) for an arriving wavefront seen in the sampling 
window. In particular, as step 7-3 correlator output analyzer 62-2A chooses a sampling 
window time index t_max at which |Y(n,t)| max occurs to be the time of arrival of the 
arriving wavefront. In addition, as step 7-4 correlator output analyzer 62-2A chooses 
5 the frequency index n_max at which |Y(n,t)| max occurs to represent the direction of 
arrival (DOA) of the arriving wavefront. The frequency index corresponds to a 
direction of arrival (e.g., 0). An amplitude for the arriving wavefront is determined as 
correlator output analyzer 62-2A divides |Y(n,t)| max by the number of antennas 
comprising the antenna array (as step 7-5). 

10 [00096] Expression 1 and the steps of Fig. 7 represent a generic non-parametric FFT 
calculation. In a CDMA-specific situation which utilizes a coding generator (such as 
coding generator 30 of Fig. 1), a comparable FFT calculation can be made using a 
refinement of Expression 1 which appears as Expression 2. 

Y(n,t) = 2 Cj *FFT(n,X(:,t)), j = 1 ,K Expression 2 

15 Expression 2 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 

[00097] As a result of operation of joint searcher and channel estimator 24-2A, an 
accurate channel estimate can be provided to the detector as a spatial signature. The 
spatial signature includes the time of arrival (TOA), as well as the direction of arrival 
20 (DOA) and amplitude. As explained below, the channel coefficient (CC) for each 
wavefront is derived from the direction of arrival (DOA) and amplitude. The time of 
arrival (TOA) and channel coefficient (CC) are applied to the detector as represented by 
lines 66-2A and 68-2A, respectively, in Fig. 2 A, 

[00098] As mentioned above, the channel coefficient (CC) for each wavefront is 
25 derived from the direction of arrival (DOA) and amplitude. Recall that at step 7-4 

correlator output analyzer 62-2A chose the frequency index n_max at which |Y(n,t)| max 
occurs to represent the direction of arrival (DOA) of the arriving wavefront, with the 
chosen frequency index corresponding to a direction of arrival (e.g., 0). The channel 
impulse response vector (i.e., array propagation vector) x is therefore generated by 
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detector interface 64-2A in accordance with Expression 3 (for identical isotropic 
antenna elements). 

[00099] x =[l,e akd * sine) , e akd * 2sin6) , ... e ^ d * (K - 1)sin9) ]*C Expression 3 

[000100] In Expression 3, j is the conventional imaginary notation; k = 2*7tX; d is a 
5 spacing distance between elements of the antenna array; X is the wavelength of the 
received/transmitted electro-magnetic signal: (f*X =c) and, K is the antenna element 
index (illustrated as antenna numbers Al, A2, A3, A4 in FIG. 9 A, for example) In 
Expression 3, C is a complex constant in which |C| = |FFT_max [/number of antennas; 
the argument of C, i.e., arg(C) = arg(FFT_max), wherein |FFT_max is the FFT value 
10 computed at step 7-1 of Fig. 7. 

[000101] In the foregoing description it is the role of channel estimate (CE) 
generator 60-2 A, and particularly detector interface 64-2 A, to generate both a time of 
arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 
from the direction of arrival, e.g., as above described in conjunction with Expression 3. 

15 In an alternate implementation of this and other embodiments described herein, the 
detector itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of 
arrival (TOA) and direction of arrival (DO A) for each arriving wavefront, may have the 
intelligence to compute the channel coefficient for each wavefront from the 
corresponding direction of arrival (DO A) information. In such case, the time of arrival 

20 and direction of arrival are output by detector interface 64 to the detector. 

[000102] Thus, considering the aspects above discussed, the joint searcher and 
channel estimator 24-2A looks in a discrete number of possible directions of arrival, 
and picks the direction of arrival with the highest correlation (highest absolute value). 
A comparative operational evaluation was preformed to illustrate the efficacy of a joint 

25 searcher and channel estimator such as joint searcher and channel estimator of 24-2A of 
Fig. 2A. A first scenario of the comparative operational evaluation involved a 
conventional searcher functioned essentially in prior art fashion for a sampling window. 
In so doing, with respect to each antenna for the sampling window the conventional 
search merely picked the time (e.g., chip) which had the greatest absolute value. In 

30 other words, the signals from each antenna were processed separately. A second 

scenario of the comparative operational evaluation was performed in the manner above 
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described with respect to the joint searcher and channel estimator 24-2A of Fig. 2 and 
Expression 1 . The same signal was applied in both scenarios to an antenna array 
having eight antenna elements. The length of the sampling window for both scenarios 
was twenty chips, and a coding sequence of { 1 } was utilized (e.g., only one of the chips 
5 contained the signal, the remainder of the chips contained complex white noise). 

[000103] Fig. 8 A illustrates the first scenario which utilized the conventional 
searcher. In contrast, Fig. 8B illustrates the spatial joint searcher and channel estimator 
24-2A of Fig. 2A utilized for the second scenario. The superiority of the second 
scenario (and thus the spatial joint searcher and channel estimator) is evident by a 
10 comparison of Fig. 8A and Fig. 8B, due to the higher SNR for the signal of interest in 
Fig. 8B . In the second scenario, it is much easier to pick out the tone or value for the 
arriving wavefront. For the second scenario, Fig. 8C(1) shows the absolute value of the 
complex channel impulse response taps; Fig. 8C(2) shows phase errors of the complex 
channel impulse response taps; and Fig. 8C(3) shows the detected time of arrival. 

15 [000104] Whereas the joint searcher and channel estimator of Fig. 2A includes a non- 
parametric type matrix analyzer, e.g., a correlator (e.g., a filter which performs a Fast 
Fourier Transform (FFT) calculation), in other example embodiments the matrix 
analyzer of the joint searcher and channel estimator implements parametric techniques. 
As does the Fig. 2A embodiment, the spatial joint searcher and channel estimator 24-2B 

20 of Fig. 2B (which uses a parametric technique) is shown along with its associated 

example antenna array 22-2B. Again by way of example, antenna array 22-2B includes 
four antenna elements 22-2B-1 through 22-2B-4. The signals obtained from the 
antenna elements are each applied to joint searcher and channel estimator 24-2B, as 
well as to a detector (not illustrated in Fig. 2B). 

25 [000105] Similar to the earlier described embodiment, joint searcher and channel 

estimator 24-2B can comprise an antenna signal matrix handling unit 40-2B, which in 
turn comprises antenna signal matrix generator 42-2B and antenna signal matrix 
memory 44-2B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-2B 

30 can also be conceptualized as matrix 80, and as such has a sampling window time 

index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 6, and is now also discussed with reference to Fig. 9A for sake of expounding the 
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joint searcher and channel estimator24-2B of Fig. 2B. 

[000106] The joint searcher and channel estimator 24-2B further comprises a matrix 
analyzer, e.g., parametric estimator 51-2B, which utilizes a parametric technique. In 
addition, in similar manner as the preceding embodiment, joint searcher and channel 
5 estimator 24-2B comprises a channel estimate generator 60-2B which has parametric 
estimation output vector analyzer 62-2B and a demodulator interface 64-2B. Basic 
steps performed by parametric estimator 51 -2B and parametric estimation output vector 
analyzer 62-2B of the joint searcher and channel estimator 24-2B of Fig. 2B are 
illustrated in Fig. 10. 

10 [000107] For each sampling window time index of the antenna signal matrix 80, as step 
10-1 the parametric estimator 51 -2B estimates, e.g. two parameters at each time instant: 
a spatial frequency parameter parameter and a spatial amplitude parameter. The spatial 
frequency parameter estimates the frequency the incident waves creates when arriving 
at the ULA. The spatial amplitude parameter estimates the amplitude of this frequency. 

15 The spatial frequency parameter and spatial amplitude parameter are considered to be a 
parameter pair and in Fig. 9B, they are illustrated as one parameter per sample along 
the sampling time index. The parameters can be calculated by an appropriate strategy 
or goal criteria, e.g., by a minimum mean square error technique (MMSE). 

[000108] As step 10-2, parametric estimation output vector analyzer 62-2B finds certain 
20 "qualifying" values in parametric estimation output vector, i.e. high or maximum values 
of the spatial amplitude parameter. The qualifying values can be, for example, values 
whose absolute values are sufficiently high or are a maximum. Each qualifying value 
of parametric estimation output vector 90 can correspond to an arriving wavefront for 
the sampling window. 

25 [000109] For each qualifying value, as step 10-3 the parametric output estimation 
vector analyzer 62-2B chooses a time of arrival (TOA) as corresponding to the 
sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs. 

[0001 10] Similarly, for each qualifying value, as step 10-4 the analyzer 62-2B chooses 
30 a direction of arrival (DOA) as the spatial frequency parameter value at the time of 
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arrival, decided in 10-3. 

[0001 11] As step 10-5, the parametric estimation output vector analyzer 62-2B 
determines the amplitude as the value of the spatial amplitude value divided by the 
number of antenna elements in the array. 

5 [0001 12] The joint searcher and channel estimator 24-2B thus looks for an 
optimum direction, and prepares a channel estimate which can be provided to the 
detector as a spatial signature. The spatial signature includes the direction of arrival 
(DOA) and amplitude. The channel coefficient (CC) for each wavefront is derived 
from the direction of arrival (DOA) and amplitude in the manner explained above with 
10 reference to Expression 3. The time of arrival (TOA) and channel coefficient (CC) are 
applied to the detector as represented by lines 66-2B and 68-2B, respectively, in Fig. 
2B. 

[0001 13] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 

15 with reference to the parametric estimation output vector 90 of Fig. 9B, the parametric 
estimation output vector analyzer 62-2B may see other high numbers and for each of 
those high numbers which qualify, an arriving wavefront may be ascertained. For 
example, if there were two high numbers, then the channel impulse response may 
reflect two arriving wavefronts. For each of the two arriving wavefronts the joint 

20 searcher and channel estimator would pick out both a time of arrival (TOA) and 

direction of arrival (DOA), as well as amplitude, which are mapped to two different 
channel coefficients, with these two different channel coefficients forming part of the 
channel estimate. 

[0001 14] Fig. 4 showed that a wavefront individually reached each of four example 
25 antenna elements of an antenna array, providing a different antenna output (complex 
baseband signal) for each antenna element. For example, the output of antenna 
element 22-1 with the complex vector a M (and phase 0m); the output of antenna 
element 22-2 is the complex vector ai_ 2 (and phase ©1.2), and so forth. The linear 
combination, of the complex antenna baseband signal and the antenna weight vectors 
30 Wi have the effect of a summation, or coherent combination in the time and space 
domain, shown as summation function 100 in Fig. 12. 
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[0001 15] The coherent combination facilitated by the antenna weight vectors Wj is 
illustrated in Fig. 11. In the example case of the four antenna elements shown in Fig. 
12, the effect of weight belonging to antenna index 2, here denoted asW 2 is to rotate the 
output of antenna element 22-2 so that its phase Q U2 lines up with the phase 0 M of the 
5 output of antenna element 22- 1 , in the manner shown in Fig. 1 1 . Similarly, the effect of 
weight W 3 is to rotate the output of antenna element 22-3 so that its phase 0i. 3 lines up 
with the phase 0 M of the output of antenna element 22-1. The effect of weight W 4 is 
to rotate the output of antenna element 22-4 so that its phase 0!_ 4 lines up with the phase 
0i_i of the output of antenna element 22-1. For simplicity, Fig. 11 ignores noise 
10 considerations, which tend to make the resultant vector less than straight. Note that in 
the preciding paragraph, the weight vectors are denoted with Wi, where i denotes the 
antenna index of the weight vector W, which is denoted without index. 

[0001 16] In the spatial joint searcher and channel estimators, the SINR for finding 
channel taps (peaks) should be proportional to the number of antenna elements 
15 comprising the array. The operation of the spatial joint searcher and channel estimators 
can be adapted to take into consideration channel variations over time, e.g., spatial 
variations in the environment (e.g., in the sending and receiving antennas). 

[0001 17] The non-parametric FFT-type correlator and the parametric techniques 
illustrated above, e.g., by Fig. 2A and Fig. 2B, respectively, are only two example 
20 techniques for finding the values or "tones" in antenna signal matrix 80 which are 

associated with arriving wavefronts. Other parametric approaches are described by or 
understood from Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis, 
ISBN-013-258419-0, Prentice Hall, which is incorporated by reference in its entirety, 
particularly Chapter 4 thereof. 

25 [0001 18] The spatial joint searcher and channel estimator and techniques of 

operation thereof as described above are suitable for any receiver unit which has plural 
receiving antennas. Thus, the spatial joint searcher and channel estimator is particularly 
well suited for, but not limited to, a base station which has plural antennas. Also 
encompasses are mobile terminals which have plural antennas. 
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[0001 19] TEMPORAL JOINT SEARCHERS/ESTIMATORS 



[000120] In other embodiments, the joint searcher and channel estimator 
contemporaneously processes the signals received at an antenna element from plural, 
successive sets of pilot data (each set of pilot data being received in its own sampling 
5 window) in order to determine both time of arrival and the channel coefficient. In so 
doing, the joint searcher and channel estimator takes into consideration a doppler shift 
or frequency shift (the terms "doppler shift" and "frequency shift" being used 
interchangeably in conjunction with the description of the temporal joint searcher and 
channel estimator). The frequency shift is primarily attributable to a doppler shift, but 
10 can also include a frequency shift in the transmitter and receiver oscillators. For 

simplification, such frequency shifts are hereinafter referred to as "doppler shifts" or 
"doppler frequency shifts". 

[000121] The doppler shift can be occasioned by movement such as relative 
movement of one of the transmitter and the receiver (for example, by movement of a 
15 mobile terminal), or movement of a signal path-affecting object or structure in the 
surroundings (which can cause a doppler shift even for a fixed transmitter and fixed 
receiver). 

[000122] In providing the channel estimate, the joint searcher and channel 
estimator essentially concurrently considers plural signals (e.g., plural sets of pilot data) 
20 received by the antenna element. The joint searcher and channel estimator applies the 
channel coefficient and the time of arrival to a detector which provides, e.g., a symbol 
estimate. 

[000123] In these embodiments, the joint searcher and channel estimator is 
essentially a two dimensional unit, with a second dimension being a temporal 

25 dimension imparted by the time intervals at which the successive sets of pilot data 
arrive. This temporal dimension, which involves essentially simultaneous and 
concurrent processing together of signals received at the antenna element from each of 
the plural sets of pilot data, bestows on these embodiments of the joint searcher and 
channel estimator the distinction of being a "temporal" joint searcher and channel 

30 estimator. 
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[000124] The temporal joint searcher and channel estimator may take differing 
embodiments and have differing implementations. In one example, illustrative 
embodiment the temporal joint searcher and channel estimator includes a non- 
parametric type correlator (e.g., a correlator which performs a Fast Fourier Transform 
5 (FFT) calculation). In another example, illustrative embodiment the temporal joint 
searcher and channel estimator utilizes a parametric approach. 

[000125] Fig. 13A illustrates one example embodiment of a spatial joint searcher and 
channel estimator 24-13A which uses a non-parametric technique for determining time 
of arrival and channel estimate, as well as an associated example antenna array 22- 13 A. 
10 In the example of Fig. 13 A, the antenna array 22-1 3 A is shown as having one antenna 
element 22-13A-1. As explained hereinafter, complex baseband signals obtained from 
the same antenna element (e.g., antenna element 22-13A-1) upon receipt of each of 
successive sets of pilot data (as hereinafter described) are each applied to joint searcher 
and channel estimator 24-13A, as well as to a detector (not illustrated in Fig. 13A). 

15 [000126] The joint searcher and channel estimator 24-1 3A comprises an antenna signal 
matrix handling unit 40- 13 A. In one particular example manifestation, antenna signal 
matrix handling unit 40-1 3 A includes antenna signal matrix generator 42-1 3 A and 
antenna signal matrix memory 44- 13 A. A matrix analyzer, which for the non- 
parametric technique of Fig. 2 A can be correlator 50- 13 A, operates on complex values 

20 stored in antenna signal matrix memory 44- 13 A. The correlator 50-1 3 A preferably 
comprises a filter. The correlator 50-13A generates certain output values, which may 
be stored, e.g., in correlator output value memory 52- 13 A. The joint searcher and 
channel estimator 24-13A further comprises a channel estimate (CE) generator 60- 13 A. 
In the illustrated example embodiment, the channel estimate (CE) generator 60-13A 

25 comprises a correlator output analyzer 62-13A and a detector interface 64- 13 A. The 
detector interface 64-13A generates, for each wavefront, a channel estimate which 
includes both a time of arrival (TO A) and a channel coefficient (CC). In Fig. 13 A, the 
time of arrival and channel coefficient output by detector interface 64-13A are applied 
to the detector on lines 66-13A and 68- 13 A, respectively. 

30 [000127] As shown in Fig. 14, the temporal joint searcher and channel estimators such 
as joint searcher and channel estimator 24-13A of Fig. 13A watch the channel response 
from an antenna (e.g., antenna 22-13-1) for sets of pilot data which are interpersed or 
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otherwise transmitted with other data (e.g., user data). For sake of simplicity, it is 
assumed that each set of pilot data is received in a separate sampling window. Such 
need not be the case, however, as differing sets of pilot data can be received 
simultaneously if the different streams are, e.g., code multiplexed. Merely as an 
5 illustrative example, Fig. 14 shows four sets of pilot data, i.e., pilot sets Tl - T4, 

interspersed with user data and received at unique global times (as indicated by the "T" 
axis in Fig. 4). 

[000128] Each set of pilot data is typically in a different frame from another set of pilot 
data. For example, pilot set Tl may be in frame 1; pilot set T2 may be in frame 11; 
10 pilot set T3 may be in frame 21; etc. "Frame transmission interval" refers to the time 
between two successive frames which contain pilot data. The time between two 
successive frames which contain pilot data is typically specified by a standard or other 
specification. 

[000129] Fig. 14 thus reflects the typical periodic transmission of pilot data by the 
15 transmitter source, and also the expected receipt of repetitions of the pilot data at the 
receiver at successive intervals. In view of factors such as relative motion of the 
transmitter and receiver, the successive intervals between differing sets of pilot data are 
not necessarily constant. 

[000130] As further shown in Fig. 14, an antenna matrix handling unit (such as antenna 
20 matrix handling unit 40-13A of the Fig. 13A embodiment) samples the signals received 
by the antenna element for each of the successive sets of pilot data, i.e., for pilot sets Tl 
- T4. Using the sampled signals, antenna signal matrix generator 42-13A generates an 
antenna signal matrix such as antenna signal matrix 110 illustrated in Fig. 14. The 
antenna signal matrix 110 may be stored in any convenient fashion, such as antenna 
25 matrix memory 44- 1 3 A. 

[000131] The antenna signal matrix 1 10 is a two dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 1 10 as a 
function of two different indexes. For the antenna signal matrix 110 shown in Fig. 14, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 14. 
30 For embodiments which utilize spreading codes or similar codes, the first index may be, 
for example, a chip index. Thus, the sampling window time index points to a time in 
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the sampling window relative to a start of the respective sampling window. In the 
antenna signal matrix 1 10 of Fig. 14, a second index, shown along the Y axis, is a pilot 
set index (which serves as a dimensional differentiation index). The pilot set index 
indicates which one of the sets of pilot data the sample was obtained. In other words, a 
5 pilot set index = Tl indicates that the sample was obtained from pilot set Tl; a pilot set 
index = T2 indicates that the sample was obtained from pilot set T2; and so on as 
depicted by the arrows which connect the matrix 110 with the received signal depiction 
with its illustrative successive sets of pilot data. As can be seen, the pilot set index 
points to a different row of the antenna signal matrix 110, each row being associated 
10 with a different set of pilot data. 

[000132] Fig. 14 shows four rows in antenna signal matrix 1 10 for consistency with the 
illustrated example wherein the antenna signal matrix encompasses four successive sets 
of pilot data. The number of sets of pilot data subsumed in a given antenna signal 
matrix, and thus the maximum value of the pilot set index, can vary from receiver to 

15 receiver, so that the present example's choice of four sets of pilot data is only 

illustrative for sake of example. In general, the choice of the number of sets of pilot 
data to be apprehended simultaneously by a temporal joint searcher and channel 
estimator depends on how quickly the doppler is expected to change. The number of 
taps/incident waves depends on the multipath. In other words in an open space we have 

20 one direct path and thus only one channel/tap coefficient in the channel impulse 
response. 

[000133] The antenna signal matrix 1 10 is also conceptualized as storing 
"dimensionally differentiated" signals acquired from a single antenna element of the 
antenna array. For the temporal joint searcher and channel estimator, wherein the 
25 antenna structure comprises an antenna which provides signals for each of successive 
sets of pilot data received at separated time intervals, the signals acquired by the 
antenna are dimensionally differentiated with regard to a temporal or time dimension. 
For example, the signals acquired by the antenna are dimensionally differentiated by 
being acquired in differing frame transmission intervals. 

30 [000134] For sake of simplicity, the complex values stored in antenna signal matrix 

110, including the complex values obtained from the antennas, are not illustrated in Fig. 
14. Such complex values would be illustrated in a third dimension, e.g., out of the 
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plane of Fig. 14. The antenna signal matrix 110 includes both complex white noise and 
(for the sake of the present illustration) a complex sample for at least one wavefront 
(planar or other known shape). The wavefronts have known phase (temporal, non- 
coherent detection), and are modulated code sequences. 

5 [000135] The complex values stored for each column of the antenna signal matrix 
1 10 of Fig. 14 can be conceptualized as a dimensional receptivity vector. That is, a 
dimensional receptivity vector is formed with the complex values taken with respect to 
a same single sampling window time index for each of the sets of pilot signals included 
in the sampling window (e.g., for sets Tl - T4 in Fig. 14). Each element taken from a 

10 unique row of antenna signal matrix 110 has a different phase in the manner of the 
differing 0 values illustrated in Fig. 5. As received by the differing antenna elements, 
for the temporal joint searcher and channel estimator the change in phase over time is 
the doppler frequency for the dimensional receptivity vector. The phase rotation speed, 
or frequency, of the dimensional receptivity vector, for the sampling window time 

15 instance can be interpreted as the doppler shift (DS). Thus, each dimensional 

receptivity vector corresponds to a separate doppler shift frequency. There are plural 
possible frequencies for the dimensional receptivity vector, each of the plural possible 
frequencies corresponding to a different possible doppler shift for a wavefront. For the 
non-parametric techniques herein employed, the plural possible frequencies can be a 

20 continuous range of frequencies. For sake of differentiating the plural possible 

frequencies, the plural possible frequencies are each represented by a frequency index. 

[000136] For the temporal joint searcher and channel estimator, the channel 
estimate comprises, as mentioned before, a time of arrival (TO A) and doppler shift for 
each arriving wavefront in the sampling window (e.g., a channel coefficient mapped to 

25 a doppler shift). Therefore, the channel estimate may comprise a set (of one or more) 
pairs of data, each pair including a time of arrival (TO A) and a channel coefficient. 
The task for the temporal joint searcher and channel estimator is thus to locate a value 
or "tone" in antenna signal matrix 110 that best corresponds to an arriving wavefront, 
e.g., to locate a value or tone for each arriving wavefront in the sampling window. This 

30 task of locating a value or "tone" in antenna signal matrix 110 that best corresponds to 
an arriving wavefront can be accomplished by various techniques, including both 
parametric and non-parametric techniques. A Fast Fourier Transform (FFT) technique 
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as discussed below is just one representative and illustrative example non-parametric 
type of correlator which can be utilized. 

[000137] Fig. 15 depicts example basic steps performed by an example correlator 
50-1 3 A and correlator output analyzer 62-13A in conjunction with the Fast Fourier 
5 Transform (FFT) calculation. As step 15-1, the correlator 50-13A of Fig. 13A 
calculates Expression 5. 

Y(n,t) = FFT(n,X(n,t)) Expression 5 

wherein t is the sampling window time index; X(n,t) is the complex antenna matrix; 
and, n is the doppler frequency index. Each FFT calculation is thus a one dimensional 
10 FFT calculation on the baseband signal, and corresponds to a specific doppler shift 
frequency. 

[000138] The output of correlator 50- 13 A, i.e., the Y(n,t) values computed using 
Expression 1, are stored as correlator output values. The correlator output values can 
be stored, for example, in the correlator output value memory 52-13A of Fig. 13 A. 

[000139] The correlator output analyzer 62-13A of channel estimate (CE) generator 
60-13A searches the correlator output values and (as step 15-2) determines therefrom a 
maximum absolute value |Y(n,t)| max . This maximum absolute value |Y(n,t)| max is 
utilized by correlator output analyzer 62-13A to determine both the doppler shift (DS) 
and time of arrival (TO A) for an arriving wavefront. In particular, as step 15-3 the 
correlator output analyzer 62-13A chooses a sampling window time index t_max at 
which |Y(n,t)| 

max occurs to be the time of arrival of the arriving wavefront. In addition, 
as step 15-4 the correlator output analyzer 62-13A chooses the doppler index n_max at 
which |Y(n,t)| max occurs to determine the doppler shift (DS) of the arriving wavefront. 
An amplitude for the arriving wavefront is determined as correlator output analyzer 62- 
13A divides |Y(n,t)| max by the number of sets of pilot data comprising the antenna signal 
matrix (as step 15-5). 

[000140] Expression 5 and the steps of Fig. 15 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
30 of Expression 5 such as that which appears as Expression 2, previously discussed, but 
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applied for the temporal joint searcher and channel estimator rather than for the spatial 
joint searcher and channel estimator. 

[000141] As a result of operation of joint searcher and channel estimator 24- 13 A, 
an accurate channel estimate can be provided to the detector as a temporal signature. 
5 For each wavefront, the temporal signature includes the time of arrival (TOA) mapped 
to a doppler (frequency) shift. As explained below, the channel coefficient (CC) for 
each time of arrival and wavefront is derived from the doppler frequency shift. The 
time of arrival (TOA) and channel coefficient (CC) are applied to the detector as 
represented by lines 66-13A and 68- 13 A, respectively, in Fig. 13 A. 

10 [000142] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the doppler frequency shift (DS). Recall that at step 15-4 the correlator 
output analyzer 62-2B chose the frequency index n_max at which |Y(n,t)| max occurs to 
represent the doppler shift frequency (DSF) of the arriving wavefront, with the chosen 
frequency index corresponding to a doppler shift (e.g., 0', i.e., the derivative of 0). The 

15 channel impulse response vector (i.e., array propagation vector) x is therefore generated 
by detector interface 64-2B in accordance with Expression 6. 

[000143] C[j 2nn + H) , e* 27Cfr2 + H) , e*^ 3 + H) , . . .J 27 ™ + H) ] Expression 6 

[000144] In Expression 6, C is the amplitude of the wavefront, f is the frequency of 
the signal (including doppler shift); T is the period time between two pilot 
20 symbols/sequences (which are assumed to be periodical, in analogy to the uniform 

array of the spatial embodiment), and H is a complex value of the signal at the first pilot 
symbol/sequence, H being arg(FFT max). For sake of simplicity, noise has been 
excluded from Expression 6, and C is assumed to be constant within the time TN. 

[000145] In the foregoing description it is the role of channel estimate (CE) 
25 generator 60-2A, and particularly detector interface 64-2A, to generate both a time of 
arrival (TOA) and a channel coefficient (CC), the channel coefficient being derived 
from the doppler shift, e.g., as above described in conjunction with Expression 6. In an 
alternate implementation of this and other embodiments described herein, the detector 
itself (such as detector 26 illustrated in Fig. 1), upon receiving the time of arrival 
30 (TOA) and doppler shift (DS) for each arriving wavefront, may have the intelligence to 



34 

compute the channel coefficient for each wavefront from the corresponding direction of 
arrival (DOA) information. In such case, the time of arrival and direction of arrival are 
output by detector interface 64-1 3 A to the detector. 

[000146] Thus, the joint searcher and channel estimator 24-1 3A looks at a discrete 
5 number of possible doppler frequency shifts, and picks the doppler frequency with the 
highest correlation (highest absolute value). 

[000147] Whereas the joint searcher and channel estimator of Fig. 13A includes a non- 
parametric correlator (e.g., a filter) which performs a Fast Fourier Transform (FFT) 
calculation, in other example embodiments the temporal joint searcher and channel 
10 estimator implements parametric techniques. As does the Fig. 13A embodiment, the 
spatial joint searcher and channel estimator 24-13B of Fig. 13B is shown along with its 
associated example antenna array 22-13B comprising an antenna element 22-13B-1 
which receives the successive sets of pilot data in the manner of Fig. 14. 

[000148] Similarly to the earlier described embodiment, joint searcher and channel 
15 estimator 24-13B can comprise an antenna signal matrix handling unit 40-13B, which 
in turn comprises antenna signal matrix generator 42-13B and antenna signal matrix 
memory 44-13B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-13B 
can also be conceptualized as matrix 110, and as such has a sampling window time 
20 index. The antenna signal matrix 110 has been previously discussed in conjunction 
with Fig. 14, and is now also discussed with reference to Fig. 16A for sake of 
expounding the joint searcher and channel estimator 24-13B of Fig. 13B. 

[000149] The joint searcher and channel estimator 24-13B further comprises the 
parametric estimator 51-13B which outputs a parametric output estimation vector for 

25 storage in memory 52-13B. In addition, in similar manner as the preceding 

embodiment, joint searcher and channel estimator 24-13B comprises a channel estimate 
generator 60-13B which has parametric output estimation vector analyzer 62-13B and a 
demodulator interface 64-13B. Basic steps performed by parametric estimator 51-1 3B 
and parametric output estimation vector analyzer 62-13B of the joint searcher and 

30 channel estimator 24-13B of Fig. 13B are illustrated in Fig. 17. 
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[000150] For each sampling window time index of the antenna signal matrix 110. As 
step 17-1, the parametric estimator 51-13B estimates, e.g., two parameters at each time 
instant: a temporal frequency parameter parameter and a temporal amplitude 
parameter. The temporal frequency parameter estimates the frequency the incident 
5 waves creates when arriving at the antenna for the consecutive pilot symbols. The 

temporal amplitude parameter estimates the amplitude of this frequency. The temporal 
frequency parameter and temporal amplitude parameter are considered to be a 
parameter pair and in Fig. 16B , they are illustrated as one parameter per sample along 
the sampling time index. 

10 [000151] As step 17-2 performed by joint searcher and channel estimator 24-13B, 
analyzer 62-13B finds certain "qualifying" values in parametric output estimation 
vector 120, i.e. maximum value of temporal amplitude vector. Each qualifying value of 
parametric output estimation vector 120 can correspond to an arriving wavefront for the 
sampling window. 

15 [000152] For each qualifying value, as step 17-3 the parametric output estimation 
vector analyzer 62-13B chooses a time of arrival (TO A) as corresponding to the 
sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs 

[000153] Similarly, for each qualifying value, as step 17-4 the parametric output 
20 estimation vector analyzer 62-13B chooses a a doppler shift frequency (DS) as the 
temporal frequency parameter value at the time of arrival decided in 17-3. 

[000154] As step 17-5, the parametric estimation output vector analyzer 62-13B 
determines the amplitude as being the maximum/qualifying absolute value divided by 
the number of pilot data sets in the series. 

25 [000155] The joint searcher and channel estimator 24-13B thus looks for an 
optimum doppler (shift) frequency, and prepares a channel estimate which can be 
provided to the detector as a temporal signature. The temporal signature includes the 
time of arrival (TOA), as well as the doppler shift frequency (DSF) and amplitude. The 
channel coefficient (CC) for each time of arrival and wavefront is derived from the 

30 doppler shift (DS) in the manner described above with reference to Expression 6. The 
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time of arrival (TO A) and channel coefficient (CC) are applied to the detector as 
represented by lines 66-13B and 68-13B, respectively, in Fig. 13B. 

[000156] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 

5 with reference to the parametric output estimation vector 120 of Fig. 16B, the 

parametric output estimation vector analyzer 62-13B may see other high numbers and 
for each of those high numbers which qualify, an arriving wavefront may be 
ascertained. For example, if there were two high numbers, then the channel impulse 
response may reflect two arriving wavefronts. For each of the two arriving wavefronts 

10 the joint searcher and channel estimator would pick out both a time of arrival (TO A) 
and doppler shift frequency (DSF), as well as amplitude, which are mapped to two 
different channel coefficients, with these two different channel coefficients forming 
part of the channel estimate for the channel impulse response. 

[000157] The operation of the temporal searcher and channel estimator has been 
15 described above for one antenna element of an antenna array 22. It should be 

understood that the antenna array 22 may comprise plural antenna elements, and that 
the operations described above may be performed separately with respect to one or 
more antenna elements of the array. Moreover, as described later, principles of the 
foregoing operation may be performed in a combined manner with respect to plural 
20 antennas of the antenna array. 

[000158] The temporal joint searcher and channel estimator and techniques of 
operation thereof as described above is particularly well suited for, but not limited to, a 
receiver unit which has only one antenna element, e.g., a mobile terminal with only one 
antenna. As indicated above, however, the temporal joint searcher and channel 
25 estimation techniques can be utilized separately but in parallel by plural antennas for a 
receiver. 

[000159] Consider, for example, the situation reflected in Fig. 1 1 in which the output of 
antenna element 22-13A-1 (or 22-13B-1) for a pilot data set Tl with the complex vector 
a M (and phase 0m); the output of the same antenna element for a pilot data set T2 with 
30 the complex vector ai_ 2 (and phase 0i_ 2 ), and so forth. In this situation, the linear 

combinationof the complex antenna baseband signal and the doppler weight vectors Wj 
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also has the effect of a summation, or coherent combination in the time domain, shown 
as summation function 100 in Fig. 12. By adding these complex vectors coherently, the 
temporal joint searcher and channel estimator increases the performance of the search 
and the channel estimate. 

5 [000160] In situations in which there is no doppler shift (e.g., the mobile terminal 
stands still or moves in a radial direction relative to the base station), the doppler shift 
frequency may be zero. In such cases the pilot data of the arriving wavefront(s) have 
essentially the same complex values. The situation of no doppler shift is just one 
special case of the generic operation of the temporal joint searcher and channel 

10 estimator described above. When the mobile starts to move, a doppler shift may occur, 
the temporal joint searcher and channel estimators obtains the doppler shift frequency, 
and thereby enhances the channel estimate. The channel estimate is enhanced by 
considering the doppler shift, regardless of the magnitude of the doppler shift. 

[000161] The non-parametric FFT-type correlator and the parametric estimator 
15 techniques illustrated above, e.g., by Fig. 13A and Fig. 13B, respectively, are only two 
example techniques for finding the values or "tones" in antenna signal matrix 110. 
Other parametric approaches are described by or understood from Stocia, Petre and 
Moses, Randolph, Introduction To Spectral Analysis, ISBN-01 3-258419-0, Prentice 
Hall, which is incorporated by reference in its entirety, particularly Chapter 4 thereof. 

20 [000 1 62] SPATIAL-TEMPORAL JOINT SEARCHERS/ESTIMATORS 

[000163] In some further embodiments, which combine features from both the 
spatial and temporal embodiments discussed above, plural antenna elements of an 
antenna array provide respective plural series of signals for successive sets of pilot data. 
The joint searcher and channel estimators of these further embodiments essentially 
25 concurrently consider the plural series of signals provided by the plural antennas for 
determining both a time of arrival and channel coefficient. 

[000164] By concurrently considering the signals provided by plural antennas, the 
channel estimate takes into consideration direction of arrival in determining the time of 
arrival and channel coefficient. By concurrently considering the series of signals 
30 provided by each antenna, in which each series comprises successive sets of pilot data, 



the channel estimate further takes into consideration a frequency shift which may be a 
doppler shift (occasioned by relative movement of a transmitter and the receiver or of 
an object in a field between the transmitter and receiver). The channel estimate is 
performed by considering spatial and temporal domain jointly and concurrently. 

5 [000165] Since it processes the series of signals from plural antennas, with each 
series comprising successive sets of pilot data, the joint searcher and channel estimator 
is considered a three dimensional unit. A first dimension is with reference to a time 
index of a sampling window, i.e., a sampling window time index. A second dimension 
is a spatial dimension imparted by the spacing of the plural antennas of the array. This 

10 spatial dimension, which involves essentially simultaneous and concurrent processing 
together of signals from the plural antennas for the array in order to determine the time 
of arrival and channel coefficient, bestows on the joint searcher and channel estimator 
the distinction of being a "spatial" joint searcher and channel estimator. A third 
dimension is a temporal dimension imparted by the time interval reflected by the 

15 successive sets of pilot data. This temporal dimension, which involves essentially 
simultaneous and concurrent processing together of signals for each of the successive 
sets of pilot data in order to determine the time of arrival and channel coefficient, 
bestows on the joint searcher and channel estimator the distinction of being a 
"temporal" joint searcher and channel estimator. In view of being both a spatial and 

20 temporal joint searcher and channel estimator, the joint searcher and channel estimator 
is also referred to as a M combined ,, spatial/temporal joint searcher and channel 
estimator, or spatio/temporal joint searcher and channel estimator. 

[000166] Concurrent consideration of the plural series of signals may be either in a 
three dimensional essentially concurrent mode or a sequenced mode. The three 

25 dimensional essentially concurrent mode involves a single step determination of the 
time of arrival and channel coefficient by simultaneously considering signals from all 
antennas of the array for all of the plural series. The sequenced mode involves a two 
step determination of the time of arrival and channel coefficient. In the sequenced 
mode, a first step comprises determining a time of arrival and direction of arrival by 

30 concurrently considering the plural signals provided by the plural antennas for a first of 
the plural series. A second step of the sequenced mode comprises refining the 
estimation of the channel coefficient based on doppler shift by concurrently considering 
the elements of the plural series having the direction of arrival determined in the first 



39 



step. This procedure could also be performed the other way around: first determining 
the time of arrival and Doppler shift and then refining the channel estimate by 
concurrently considering the elements of the plural series having the Doppler shift 
determined in the first step. 

5 [000167] Fig. 18A illustrates one example embodiment of a spatio-temporal joint 
searcher and channel estimator 24- 13 A, as well as an associated example antenna array 
22- 18 A. The antenna array 22-18A includes, by way of non-limiting example, four 
antenna elements 22-18A-1 through 22-1 8 A-4. While the antenna elements 22-18A-1 
through 22- 18 A-4 are shown as forming a uniform linear array (ULA), it should be 

10 understood that antenna configurations other than a uniform linear are possible, and that 
the number of antenna elements in the antenna array may vary (e.g., the number of 
antenna elements is not limited to four). After suitable radio frequency processing, the 
signals obtained from the antenna elements are each applied as baseband signals to joint 
searcher and channel estimator 24- 18 A, as well as to a detector (not illustrated in Fig. 

15 18A). 

[000168] The joint searcher and channel estimator 24-1 8 A comprises an antenna signal 
matrix handling unit 40- 18 A. In one particular example manifestation, antenna signal 
matrix handling unit 40-18A includes antenna signal matrix generator 42-18A and 
antenna signal matrix memory 44- 18 A. A matrix analyzer, which for the non- 
20 parametric technique of Fig. 18A can be a correlator 50-1 8 A, operates on complex 
values stored in antenna signal matrix memory 44- 18 A. The correlator 50-18A 
preferably comprises a filter. The correlator 50-18A generates certain output values, 
which may be stored, e.g., in correlator output value memory 52- 18 A. The joint 
searcher and channel estimator 24-18A further comprises a channel estimate (CE) 
25 generator 60- 18 A. In the illustrated example embodiment, the channel estimate (CE) 
generator 60-18A comprises a correlator output analyzer 62-18A and a detector 
interface 64-18A. The detector interface 64-18A generates, for each wavefront, a 
channel estimate which includes both a time of arrival (TO A) and a channel coefficient 
(CC). In Fig. 18 A, the time of arrival and channel coefficient output by detector 
30 interface 64 are applied to the detector on lines 66-1 8 A and 68- 18 A, respectively. 

[000169] In the joint searcher and channel estimator 24-18A of Fig. 18 A, for each series 
of sets of pilot data (represented by pilot data sets Tl - T4), the antenna matrix handling 
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unit 40-1 8 A samples the signals from each antenna element. Using the sampled 
signals, antenna signal matrix generator 42-18A generates an antenna signal matrix 
such as antenna signal matrix 130 illustrated in Fig. 19. The antenna signal matrix 130 
may be stored in any convenient fashion, such as antenna matrix memory 44- 18 A. 

5 [000170] The antenna signal matrix 130 is a three dimensional functionally dependent 
matrix. In other words, complex samples are stored in antenna signal matrix 130 as a 
function of three different indexes. For the antenna signal matrix 130 shown in Fig. 19, 
a first index is a sampling window time index, illustrated along the X axis of Fig. 19. 
For embodiments which utilize spreading codes or similar codes, the first index may be, 

10 for example, a chip index. Thus, the sampling window time index points to a time in 
the sampling window relative to a start of the sampling window. 

[000171] In the antenna signal matrix 130 of Fig. 19, a second index, shown along the 
Y axis, is an antenna index. The antenna index points to a different row of the antenna 
signal matrix 130, each row being associated with a different antenna element in 

15 antenna array 22. Fig. 19 shows four rows in antenna signal matrix 130 for consistency 
with the previous examples of an antenna array comprising four antenna elements. It is 
reiterated, however, that the number of antennas in an antenna array, and thus the 
number of rows in antenna signal matrix 130 and the maximum value of the antenna 
index, can vary from receiver to receiver, and that the choice of four antenna is only 

20 illustrative for sake of example. 

[000172] In the antenna signal matrix 130 of Fig. 19, a third index, shown along the Z 
axis, is a pilot set index. The pilot set index indicates which one of the sets of pilot data 
the sample was obtained. In other words, a pilot set index = Tl indicates that the 
sample was obtained from pilot set Tl ; a pilot set index = T2 indicates that the sample 
25 was obtained from pilot set T2; and so on as depicted by the arrows which connect the 
matrix 110 with the received signal depiction with its illustrative successive sets of pilot 
data. As can be seen, the pilot set index points to a different plane of the antenna signal 
matrix 110, each plane being associated with a different set of pilot data. 

[000173] Fig. 19 shows four planes in antenna signal matrix 130 for consistency with 
30 the illustrated example wherein the antenna signal matrix encompasses four successive 
sets of pilot data. The number of sets of pilot data subsumed in a given antenna signal 
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matrix, and thus the maximum value of the pilot set index, can vary from receiver to 
receiver, so that the present example's choice of four sets of pilot data is only 
illustrative for sake of example. In general, the choice of the number of sets of pilot 
data to be apprehended simultaneously by a spatio-temporal spatial/temporal joint 
5 searcher and channel estimator depends on how quickly the doppler is expected to 
change. The number of taps/incident waves depends on the multipath. In other words 
in an open space we have one direct path and thus only one channel/tap coefficient in 
the channel impulse response. 

[000174] For sake of simplicity, the complex values stored in antenna signal matrix 
10 130, including the complex values obtained from the antennas, are not illustrated in Fig. 
19. Such complex values would be illustrated in a fourth dimension. 

[000175] In conjunction with the antenna signal matrix 130 of Fig. 19, and particularly 
a WCDMA case in which spacing of the antenna elements in the antenna array is not 
too far apart, the plane wavefront arriving at the antenna array can be considered to 
15 arrive in the same sampling window time index (or chip index). 

[000176] Assuming that the wavefront arrives at the antenna elements at different 
times (the time differences are small in comparison to the sampling time interval), the 
complex values stored for each column of the antenna signal matrix 130 of Fig. 19 have 
differing phase (e.g., 0) values in each row of the column. For a uniformly spaced array 

20 of antenna elements, the phase difference is essentially the same between adjacent rows 
of the same column (although noise may be a factor). But whatever the spacing, the 
rate of change of the phase with respect to time (time of travel of the approaching 
wavefront) is the phase rotation speed, or frequency, for the vector formed by the 
column, as previously explained. This per column frequency can be interpreted as a 

25 direction of arrival (DO A). There are plural possible frequencies for the columns of 

antenna signal matrix 130, with each of the plural possible frequencies corresponding to 
a possible direction of arrival (DO A) of a wavefront. The plural possible direction of 
arrival frequencies are represented by a frequency index "ni". 

[000177] In a similar manner, for each slice of antenna signal matrix 130 along the 
30 Z" direction the complex values have differing phase (e.g., 0) values. The Z-aligned 
elements of differing "Z" planes of antenna signal matrix 130 have differing phase 
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values in view of a possible doppler shift as detected by the different sets of pilot data 
as gathered over plural sets of pilot data in a series. The rate of change over time of the 
phase along the Z direction between successive sets of pilot data is a frequency which is 
associated with the doppler shift. There are plural possible frequencies for the Z slices 
5 of antenna signal matrix 130, with each of the plural possible frequencies 

corresponding to a possible doppler shift (DS) for wavefront. The plural possible 
doppler shift frequencies are represented by a frequency index "n 2 ". 

[000178] The channel estimate generator 60-18A (see Fig. 18 A) seeks to develop a 
"composite 11 channel estimate based on the complex values stored in antenna signal 

10 matrix 130. As mentioned before, since an antenna array such as antenna array 22-18A 
has plural antenna elements, there are a corresponding plurality of channels through 
which wavefronts are received, and accordingly there could also be a separate channel 
impulse response or separate channel estimate for each of the plural channels. But by 
storing the complex samples in antenna signal matrix 130 in the manner aforedescribed, 

15 and by concurrently finding the time of arrival (TO A) and channel coefficients over the 
entire antenna signal matrix 130, the channel estimate generator 60-18A provides a 
channel estimate which encompasses all channels for all antenna elements and for this 
reason is known as a "composite" channel estimate. 

[000179] The composite channel estimate comprises, as mentioned before, a time 
20 of arrival (TO A) and channel coefficient for each arriving wavefront in the sampling 
window (e.g., a channel coefficient mapped to a time of arrival (TOA)). Therefore, the 
channel estimate may comprise a set (of one or more) pairs of data, each pair including 
a time of arrival (TOA) and channel coefficient. The task for correlator 50-18A is thus 
to locate a value or "tone" in antenna signal matrix 130 that best corresponds to an 
25 arriving wavefront, e.g., to locate a value or tone for each arriving wavefront in the 
sampling window. 

[000180] The task of locating, in an antenna signal matrix such as antenna signal 
matrix 130, a value or "tone" that best corresponds to an arriving wavefront can be 
accomplished by various techniques, including both parametric and non-parametric 
30 techniques. A Fast Fourier Transform (FFT)) performed in a three dimensional 
essentially concurrent mode is discussed below in conjunction with as just one 
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representative and illustrative example of a non-parametric technique wherein 
correlator 50-18A is utilized. 

[000181] Fig. 20 depicts example basic steps performed by an example correlator 
50-1 8 A and analyzer 62-18A in conjunction with the Fast Fourier Transform (FFT) 
5 calculation. In conjunction with Fig. 20, Fig. 21 shows an antenna signal matrix; a 
doppler weight and antenna weight vector; and a non-parametric estimation output 
vector for an example embodiment of a spatio-temporal joint searcher and channel 
estimator which operates in a three dimensional essentially concurrent mode. As step 
20-1, the correlator 50-18A of Fig. 18A calculates Expression 8. 

10 Y(ni,n 2 ,t) = FFT(ni,n 2 ,X(:,:t)) Expression 8 

In Expression 8, t is the sampling window time index; X(:,:,t) is the complex antenna 
matrix (with the colon representing all antenna indexes for one sampling window 
time index); n t is the direction of arrival frequency index; and n 2 is the doppler shift 
index. Each FFT calculation is thus a two dimensional FFT calculation on the 
15 baseband signal, corresponding both to a specific direction of arrival (as depicted by the 
frequency index n x ) and to a specific doppler shift (as depicted by the frequency index 
n 2 ). 

[000182] The output of correlator 50-18A, i.e., the Y(n l9 n 2 ,t) values computed 
using Expression 8, are stored as correlator output values. The correlator output values 
20 can be stored, for example, in the correlator output value memory 52-18A of Fig. 18 A. 

[000183] The correlator output analyzer 62-18A of channel estimate (CE) generator 
60-1 8 A searches the correlator output values Y(n!,n 2 ,t) and (as step 20-2) determines 
therefrom a maximum absolute value |Y(ni,n 2 ,t)| max . This maximum absolute value 
|Y(n!,n 2 ,t)| max is utilized by correlator output analyzer 62-18A to determine both the 

25 direction of arrival (DOA) and time of arrival (TOA) for an arriving wavefront seen in 
the sampling window. In particular, as step 20-3 the correlator output analyzer 62-18A 
chooses a sampling window time index t_max at which |Y(ni,n 2 ,t)| max occurs to be the 
time of arrival of the arriving wavefront. In addition, as step 20-4 the correlator output 
analyzer 62-18A chooses the frequency index n^max at which |Y(n 1 ,n 2 ,t)| max occurs to 

30 determine the direction of arrival (DOA) of the arriving wavefront. Further, as step 20- 
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5 the correlator output analyzer 62-18A chooses the index n 2 _max at which 
|Y(n 1 ,n 2 ,t)| max occurs to determine the doppler shift of the arriving wavefront. An 
amplitude for the arriving wavefront is determined as the correlator output analyzer 62- 
18A divides |Y(n 1? n 2 ,t)| max by the product of the number of antennas comprising the 
5 antenna array and the number of sets of pilot data included in the matrix 130 (as step 
20-6). 

[000184] Expression 8 and the steps of Fig. 20 represent a generic FFT calculation. 
In a CDMA-specific situation which utilizes a coding generator (such as coding 
generator 30 of Fig. 1), a comparable FFT calculation can be made using a refinement 
10 of Expression 8 which appears as Expression 9. 

Y(n!,n 2 ,t) = 2 Cj *FFT(n 1 ,n 2 ,X(:,:,t)), j = 1,K Expression 9 

Expression 9 is understood from Expression 1, it being further mentioned that Cj is a 
coding sequence symbol value j; and K is a length of the coding sequence. 

[000185] As a result of operation of spatio-temporal joint searcher and channel 
15 estimator 24- 18 A, an accurate channel estimate can be provided to the detector as a 
spatio-temporal spatial and temporal signature. The spatial signature includes the 
direction of arrival; the temporal signature includes the doppler shift. The channel 
coefficient (CC) for each time of arrival and antenna element is derived from the 
direction of arrival (DOA)and the doppler shift. The time of arrival (TO A) and channel 
20 coefficient (CC) are applied to the detector as represented by lines 66-18A and 68- 18 A, 
respectively, in Fig. 18A. 

[000186] As mentioned above, the channel coefficient (CC) for each wavefront is 
derived from the direction of arrival (DO A) and doppler shift (DS). Recall that at step 
18-4 analyzer 62-18A chose the frequency index n x _max at which |Y(ni,n 2 ,t)| max occurs 
25 to represent the direction of arrival (DOA) of the arriving wavefront, with the chosen 
frequency index corresponding to a direction of arrival (e.g., 0). Further, analyzer 62- 
18A chose the frequency index n 2 _max at which |Y(n 1? n 2 ,t)| max occurs to represent the 
doppler shift of the arriving wavefront, with the chosen frequency index corresponding 
to a doppler shift. The channel impulse response vector (i.e., array propagation vector) 
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x is therefore generated by detector interface 64-18A in accordance with Expression 10 
(for identical isotropic antenna elements). 



e Gkd*sin9) e 0kd*2sine) e Gkd*(K-l)sinG)j*Qj. 



Q e Gkd*sinG) e (jkd*2sin9) e 0kd*(K-l)sin8)j* cl . 

5 (l,e akd * sin6) , e akd * 2sin6) , . . . e * d * (K " 1)sine) ]*CN Expression 10 

In Expression 10, CN = ei 2nfrN = H with H and other parameters being as previously 
defined. 

[000187] In the foregoing description it is the role of channel estimate (CE) 
generator 60- 18 A, and particularly detector interface 64- 18 A, to generate both a time of 

10 arrival (TO A) and a channel coefficient (CC), the channel coefficient being derived 
from the direction of arrival and doppler shift, e.g., as above described in conjunction 
with Expression 11. In an alternate implementation of this and other embodiments 
described herein, the detector itself (such as detector 26 illustrated in Fig. 1), upon 
receiving the time of arrival (TO A), the direction of arrival (DO A), and the doppler 

15 shift for each arriving wavefront, may have the intelligence to compute the channel 
coefficient for each wavefront from the corresponding direction of arrival (DOA) and 
doppler shift information. In such case, the time of arrival, direction of arrival, and 
doppler shift are output by detector interface 64 to the detector. 

[000188] The operation of the correlator 50-18A in calculating Expression 8 or 
20 Expression 9 is an example of a three dimensional essentially concurrent mode, since 
evaluation of Expression 8 (or Expression 9 for a WCDMA implementation) involves a 
single step determination of the time of arrival and channel coefficient by simultaneous 
considering signals from all antennas of the array for all of the plural series. In other 
words, in the illustrated example of the three dimensionally essentially concurrent 
25 mode, the Fast Fourier Transform (FFT) of Expression 8 or Expression 9 had three 
arguments: n { , n 2 , and X(:,:t), so that the FFT operated on all arguments essentially 
simultaneously. 
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[000189] In contrast to the three dimensional essentially concurrent mode, the 
sequenced mode involves a two step determination of the time of arrival and channel 
coefficient. In a first alternative way for implementing the sequenced mode, a first step 
comprises determining a time of arrival and direction of arrival by concurrently 
5 considering the plural signals provided by the plural antennas for a first of the plural 
series. For example, the first step of the first alternative of the sequenced mode can 
involve calculating a FFT such as that of Expression 1 (or, for WCDMA, Expression 
2). From the results of the first step or first FFT calculation, a time of arrival (TO A) 
and tentative channel coefficient are determined. Then, as a second step of the first 

10 alternative of the sequenced mode, the tentative channel coefficient is refined by taking 
into consideration a possible frequency shift (e.g., doppler shift) by further considering 
the elements of the plural series having the direction of arrival determined in the first 
step. In a second alternative way of implementing the sequenced mode, the order of the 
steps is essentially reversed: first the FFT is performed in the temporal domain to 

15 decide time of arrival and tentative channel coefficient; and secondly the tentative 
channel coefficient is refined by FFT in the spatial domain. 

[000190] The procedures of the first alternative implementation of the sequenced 
mode for the non-parametric technique are illustrated in Fig. 22A and Fig. 22B in 
conjunction with Fig. 23. Fig. 22 A and Fig. 22B as a diagrammatic view of an antenna 

20 signal matrix; an antenna weight vector; and a non-parametric estimation output vector 
for an example embodiment of a sequential spatio-temporal joint searcher and channel 
estimator. In Fig. 22A, the FFT operates on the spatial domain and calculates the FFT 
(illustrated by the FFT vector Wi) for the antenna matrix for each time interval. The 
time of arrival is choosen by picking the direction of arrival index and time index with 

25 the highest absolute value. If this index does not coincide for all time intervals, the 
index can be chosen with some method, e.g. majority decision. 

[000191] After having chosen the time of arrival index and direction of arrival 
index, these FFI -processed samples are further FFI -processed by a FFT calculation in 
the temporal domain (illustrated by FFT frequency vector Wj). Fig. 22B shows the 
30 spatially filtered samples for the identified time of arrival and direction (marked as grey 
in the figure) are filtered with the temporal vectors. After the second FFT processing, 
the channel estimate is created from the sample with the highest magnitude. Step 23-1 
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through step 23-7 of Fig. 23 also describe the procedure of the first alternative 
implementation of the sequenced mode. 

[000192] The procedures of the second alternative implementation of the sequenced 
mode for the non-parametric technique are illustrated in Fig. 24A and Fig. 24B in 
5 conjunction with Fig. 25. Fig. 24A and Fig. 24 B show an antenna signal matrix; a 
Doppler weight vector; and a non-parametric estimation output vector. In Fig. 24A the 
FFT operates on the temporal domain and calculates the FFT (illustrated by the FFT 
vector Wj) for the antenna matrix for each time interval. The time of arrival is chosen 
by picking the doppler index of and time index with the highest absolute value. If this 

10 index does not coincide for all time intervals, the index can be chosen with some 
method, e.g. majority decision. After having chosen the time of arrival index and 
doppler index, these FhT -processed samples are further FFT-processed by a FFT 
calculation in the spatial domain (illustrated by FFT frequency vector Wi). Fig. 24B 
shows the spatially filtered samples for the identified time of arrival and Doppler shift 

15 (marked as grey in the figure) are filtered with the spatial vectors. After the second 
FFT processing, the channel estimate is created from the sample with the highest 
magnitude. Step 25-1 through step 25-7 of Fig. 25 also describe the procedure of the 
second alternative implementation of the sequenced mode. 

[000193] Whereas the joint searcher and channel estimator of Fig. 18A includes a non- 
20 parametric type correlator (e.g., a filter which performs a Fast Fourier Transform (FFT) 
calculation), in other example embodiments the joint searcher and channel estimator 
implements parametric techniques. As does the Fig. 18A embodiment, the parametric 
temporal joint searcher and channel estimator 24-18B of Fig. 18B is shown along with 
its associated example antenna array 22-18B. Again by way of example, antenna array 
25 22-18B includes four antenna elements 22-18B-1 through 22-18B-4. The signals 
obtained from the antenna elements are each applied to joint searcher and channel 
estimator 24-18B, as well as to a detector (not illustrated in Fig. 18B). 

[000194] Similarly to the earlier described embodiment, joint searcher and channel 
estimator 24-18B can comprise an antenna signal matrix handling unit 40-18B, which 
30 in turn comprises antenna signal matrix generator 42-18B and antenna signal matrix 
memory 44-18B, which function much in the manner previously described. For 
example, the complex baseband values stored in antenna signal matrix memory 44-18B 
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can also be conceptualized as matrix 130, and as such has a sampling window time 
index. The antenna signal matrix 80 has been previously discussed in conjunction with 
Fig. 19. 

[000195] The joint searcher and channel estimator 24-18B further comprises a 
5 parametric estimator 51-18B which produces a parametric estimation output vector. In 
addition, in similar manner as the preceding embodiment, joint searcher and channel 
estimator 24-18B comprises a channel estimate generator 60-18B which has parametric 
output estimation vector analyzer 62-18B and a demodulator interface 64-18B 

[000196] Fig. 26 shows an antenna signal matrix and a parametric estimation output 
10 vector for an example embodiment of a spatio-temporal joint searcher and channel 
estimator. As with the non-parametric techniques, the parametric techniques can be 
implemented either in a three dimensional essentially concurrent mode or in a 
sequenced mode, with the sequence mode having two alternative implementations. 

[000197] Fig. 27 shows basic, representative steps involved in a parametric three 
15 dimensional essentially concurrent mode. Step 27-1 shows the joint searcher and 
channel estimator 24-18B producing a parametric estimation output vector. Then, as 
step 27-2, analyzer 62-18B finds the "qualifying" values in the parametric estimation 
output vector. 

[000198] For each qualifying value, as step 27-3 the parametric output estimation 
20 vector analyzer 62-18B chooses a time of arrival (TO A) as corresponding to the 

sampling window time index t for the qualifying value, e.g., the time index at which the 
maximum/qualifying absolute value of the parametric estimation output vector occurs. 

[000199] For each qualifying value, as step 27-4 the parametric output estimation 
vector analyzer 62-18B chooses a spatio-temporal frequency parameter corresponding 
25 to the spatio-temporal frequency for the maximum/qualifying absolute value of the 
parametric estimation output vector. 



[000200] As step 27-5, the parametric estimation output vector analyzer 62-13B 
determines the amplitude as the value for the spatio-temporal amplitude value for the 
time of arrival decided in step 27-2. 



49 

[000201] It should be understood from the foregoing that information indicative of 
more than one incident wavefront may be seen in a sampling window. For example, 
with reference to the parametric estimate output vector 140 of Fig. 26, the parametric 
output estimation vector analyzer 62-1 8B may see other (e.g., plural) high numbers, and 
5 for each of those high numbers which qualify, an arriving wavefront may be 
ascertained. 

[000202] The procedures of the first alternative implementation of the sequenced 
mode for the parametric technique are illustrated in Fig. 28A and Fig. 28B in 
conjunction with Fig. 29. Fig. 28A and Fig. 28B depict a parametric, sequential spatio- 

10 temporal joint searcher and channel estimator for this first alternative implementation. 
In Fig. 28A and Fig. 28B the parametric approach first operates on the spatial domain 
and calculates the spatial frequency parameters for each time instant over the time 
transmission intervals. The time of arrival is chosen by picking the spatial frequency 
amplitude value with the highest absolute value. The direction of arrival, DOA, is the 

15 value of the spatial frequency parameter. If this time of arrival does not coincide for all 
the time intervals, the time of arrival can be chosen with some method, e.g. majority 
decision. As shown in Fig. 28B, after having chosen the time of arrival index and 
direction of arrival, these samples are processed by the parametric approach applied in 
the temporal domain. After the second processing, the channel estimate is created from 

20 the temporal parameters. Step 29-1 through step 29-5 of Fig. 29 also describe the 
procedure of the first alternative implementation of the parametric sequenced mode. 

[000203] Fig. 30A and Fig. 30B show a parametric, sequential spatio-temporal 
joint searcher and channel estimator for a second alternative implementation of a 
parametric, sequential spatio-temporal joint searcher and channel estimator. In Fig. 

25 30A and Fig. 30B the paremetric approach first operates on the temporal domain and 
calculates the temporal frequency parameters for each time instant over the time 
transmission intervals. The time of arrival is chosen by picking the temporal frequency 
amplitude value with the highest absolute value. The Doppler shift frequency, DSF is 
the value of the temporal frequency parameter. If this time of arrival does not coincide 

30 for all the time intervals, the time of arrival can be chosen with some method, e.g. 
majority decision. As shown in Fig. 30B, after have chosen the time of arrival index 
and DSF, these samples are processed by parametric approach applied in the spatial 
domain. After the second processing, the channel estimate is created from the spatial 
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parameters. Step 31-1 through step 31-7 of Fig. 31 also describe the procedure of the 
second alternative implementation of the parametric sequenced mode. 

[000204] The non-parametric FFT-type correlator and the parametric linear 
combination logic techniques illustrated above are only two example techniques for 
5 finding the values or "tones" in antenna signal matrix 130 which are associated with 
arriving wavefronts. Other parametric approaches are described by or understood from 
Stocia, Petre and Moses, Randolph, Introduction To Spectral Analysis, ISBN-013- 
258419-0, Prentice Hall, which is incorporated by reference in its entirety, particularly 
Chapter 4 thereof. 

10 [000205] The spatio-temporal joint searcher and channel estimator and techniques 
of operation thereof as described above are suitable for any receiver unit which has 
plural receiving antennas. Thus, the spatial joint searcher and channel estimator is 
particularly well suited for, but not limited to, a base station which has plural antennas. 
The spatio-temporal joint searcher and channel estimator and techniques of operation 

15 thereof also encompasses mobile terminals which have plural antennas. 

[000206] The joint searcher and channel estimators thus employ a multi- 
dimensional and optimum detection and estimation approach. The multi-dimensional 
joint searcher and channel estimators typified by those described herein have better 
20 performance than traditional one dimensional searchers. The multi-dimensional joint 
searcher and channel estimators have a greater SNIR for detecting time of arrival, 
which increases the probability of that the correct time of arrival will be ascertained. 
This, in turn, leads to a better channel estimate. 

[000207] In terms of implementation, the blocks, units, and functionalities of the 
25 differing embodiments of the joint searcher and channel estimator herein described can 
various forms. For example, those skilled in the art will appreciate that one or more of 
the functionalities of the joint searcher and channel estimator can be implemented using 
individual hardware circuits, using software functioning in conjunction with a suitably 
programmed digital microprocessor or general purpose computer, using an application 
30 specific integrated circuit (ASIC), and/or using one or more digital signal processors 
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(DSPs). Furthermore, the functionalities of the joint searcher and channel estimator 
need not be delineated specifically in the manners illustrated, it being understood (for 
example) that the functionalities can be distributed, combined, subdivided, or otherwise 
rearranged for accomplishing essentially the same results. 

[000208] Use and operation of the joint searcher and channel estimators is not 
confined to WCDMA transmission, although in some instances WCDMA has been 
described above as an example environment of implementation. The principles, 
techniques, methods, and apparatus described herein can be adapted or augmented for 
compatibility with various types of networks, not only WCDMA, but other networks as 
well (such as GSM, for example). 

[000209] In the foregoing, it will be appreciated that other aspects of wireless 
receiver structure and operation which are tangential to matters described above have 
been omitted for clarity. Such aspects, well understood by persons skilled in the art, 
include without limitation pulse shaping, sampling frequency, time jitter, time 
15 alignment, demodulation, inter symbol interference (ISI), and co channel interference 
(CCI). 

[000210] While the invention has been described in connection with what is 
presently considered to be the most practical and preferred embodiment, it is to be 
understood that the invention is not to be limited to the disclosed embodiment, but on 
20 the contrary, is intended to cover various modifications and equivalent arrangements 
included within the spirit and scope of the appended claims. 
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